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SUMMARY 


This  report  presents  the  results  of  a  research  program  to 
obtain  simultaneously  a  set  of  blade  strain  and  air-load  data  for  a 
blade  of  a  stoppable  rotor  configuration  arbitrarily  oriented  with 
respect  to  the  free  stream.  The  research  program  sumntarizrd 
herein  obtained  data  for  such  a  blade  configuration  over  a  wide 
range  of  shaft  tilt  angles,  blade  pitch  angles, and  azimuth  angles 
during  a  wind-tunnel  test  program  conducted  in  the  NASA -Ames 
40 -by -80 -foot  wind  tunnel. 

The  main  results  of  the  program  were  as  follows: 

At  azimuth  angles  of  90  and  270  degrees,  the  variations 
of  the  aerodynamic  characteristics  with  angle  of  attack 
were  generally  those  that  might  be  expected  on  the  basis 
of  sectional  data. 

The  effect  of  blade  beamwisc  flexibility  is  very  signifi¬ 
cant  as  regards  the  variation  of  the  aerodynamic  char¬ 
acteristics  with  azimuth  angle. 

In  the  azimuth  angle  range  of  2255  f  &.  3  30  degrees,  an 
instability  was  encountered  for  both  positive  and  negative 
blade  pitch  angles. 

The  contribution  of  the  flexible  blade  bending  slope  to 
the  effective  angle  of  attack  was  significant  in  the  estab¬ 
lishment  of  the  instability  boundary. 

A  linearized  lifting -surface  theory  for  estimating  aero¬ 
dynamic  forces  on  the  stopped  rotor  blade  was  formu¬ 
lated.  Limited  comparisons  were  made  with  test  results 
and  with  other  theoretical  solutions. 


FOREWORD 


This  report  on  the  work  that  was  accomplished  during  the 
contract  effort  to  investigate  the  aerodynamic  and  aeroelastic  char¬ 
acteristics  of  a  full-scale  rotor  blade  stopped  in  flight  has  been 
prepared  in  two  parts.  This  report  presents  a  detailed  discussion 
of  all  aspects  of  the  program  and  the  analysis  of  the  results  that 
were  obtained.  All  of  the  basic  data  that  were  collected  during  the 
test  phase  of  the  program  are  on  file  at  the  U.  S.  Army  Aviation 
Materiel  Laboratories  (USAAVLABS). 

The  research  program  was  conducted  by  the  Cornell 
Aeronautical  Laboratory,  Inc.  (CAL)  under  USAAVLABS  Contract 
DA  44-177-AMC-366(T)  (Task  1 F162204A1  3902)  and  was  carried 
out  under  the  technical  cognizance  of  Mr.  William  E.  Nettles  of 
USAAVLABS. 

The  research  program  began  in  July  1966  and  was  completed 
in  June  1968.  Personnel  associated  with  the  research  program  included 
Richard  P.  White,  Jr.,  Eugene  Skelly,  Stephen  King,  and  Joseph  Nenni 
of  CAL.  Mr.  John  McCloud  of  the  NASA-Ames  research  staff  contrib¬ 
uted  many  valuable  suggestions  during  the  formulation  and  conduct  of 
the  tests  in  the  full-scale  tunnel  and  was  the  NASA  project  engineer 
during  the  test  phase  of  the  program. 
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INTRODUCTION 


During  the  past  decade,  the  flight  speed  obtainable  with  VTOL 
aircraft  has  more  than  doubled.  The  reason  for  this  rapid  increase  in 
flight  speed  has  been  twofold:  the  rapid  advance  in  VTOL  technology 
and  the  introduction  of  the  lightweight  free -turbine  engine.  The  latest 
advances  in  flight  have  been  made  with  compound -type  helicopters  in 
which  the  rotor  is  slowed  and  unloaded  by  means  of  a  lifting  surface. 

It  is  reasonable  to  anticipate  that  the  next  development  step  would  be 
toward  a  stopped  rotor  configuration  in  which  the  rotor  now  being  car¬ 
ried  on  a  compound  vehicle  would  be  completely  unloaded,  decelerated, 
stopped,  locked  in  position,  and  trailed  or  stowed. 

Stoppable  rotors  will  necessarily  traverse  a  peculiar  operating 
range.  During  the  last  revolution  of  the  rotor,  before  coming  to  a  com¬ 
plete  stop,  the  blades  will  effectively  vary  in  behavior  between  that  of  a 
high-aspect-ratio  wing  at  azimuth  positions  of  90  degrees  and  270  degrees 
to  a  very  low-aspect-ratio  wing  at  0  degrees  and  180  degrees.  The  aero¬ 
dynamic  load  amplitude  and  distribution  on  a  blade  would  be  expected  to 
undergo  large  variations  during  this  transition  from  "high-  to  low-aspect- 
ratio"  operation  and  strong  nonlinear  effects  can  be  anticipated  as  well  as 
the  possibility  of  aeroelastic  problems  such  as  static  divergence  and  flutter. 

There  is  no  verified  theory  for  predicting  the  aerodynamic  loads 
needed  to  study  the  blade  deformations  and  stability  of  rotors  stopped  at 
any  arbitrary  azimuth  position  in  flight.  If  stoppable  rotor  configurations 
are,  in  fact,  going  to  be  developed  on  a  consistent  and  rational  basis,  then 
there  is  a  need  for  a  proven  method  of  predicting  the  loads  developed  by 
such  a  rotor  system  in  flight.  The  primary  purpose  of  the  research  effort 
reported  herein,  therefore,  was  to  obtain  a  set  of  aerodynamic  loading 
data  that  was  consistent  with  the  resultant  blade  deformations  for  a  range 
of  pertinent  aircraft  flight  and  control  parameters  for  use  in  checking  the 
validity  of  any  theoretical  prediction  method  that  is  developed.  A  secon¬ 
dary  purpose  of  the  research  effort  was  to  develop  a  linearized  aero¬ 
dynamic  prediction  method  which  includes  the  effects  of  blade  deformation 
in  order  to  determine,  through  comparison  with  experimental  results,  the 
degree  to  which  such  a  linearized  theory  could  predict  the  measured  results. 

This  report  describes  the  test  equipment,  its  operation  and  cali¬ 
bration,  and  the  development  of  the  theoretical  prediction  method;  it 
presents  and  discusses  the  major  experimental  and  theoretical  results 
that  were  obtained  during  the  research  program.  The  detailed  pressure 
measurements,  blade  strain  data,  and  integrated  blade  loading  data  for 
the  approximately  7  50  configurations  tested  are  on  file  at  USAAVLABS. 
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DISCUSSION 


A.  EQUIPMENT  AND  TEST  APPARATUS 

The  purpose  of  the  test  program  was  to  obtain  simultaneously 
aerodynamic  and  structural  deformation  data  for  a  "typical"  full-scale 
rotor  blade  in  a  stopped  configuration.  A  special  blade -retention  struc¬ 
ture  was  constructed  that  incorporated  all  of  the  desired  test  features 
for  the  efficient  conduct  of  the  experimental  program,  and  a  set  of 
specially  instrumented  blades  that  had  been  constructed  previously 
and  tested  in  a  normal  helicopter  rotor  configuration  was  used  as  the 
basic  rotor. 

The  special  two-bladed  test  rotor  system  that  was  assembled 
for  the  program  is  shown  ready  for  test  in  the  NASA -Ames  40-by-80- 
fcot  wind  tunnel  in  Figure  1.  While  a  two-bladed  rotor  system  was  used 
for  the  tests,  only  one  blade  was  fully  instrumented  to  measure  the  chord- 
wise  pressure  distributions  at  six  spanwise  stations,  and  the  beamwise 
bending,  edgewise  bending,  and  torsional  strains  at  nine  spanwise  sta¬ 
tions.  In  addition  to  these  quantities,  the  normal  force,  beamwise 
moment,  edgewise  force,  edgewise  moment,  and  pitching  moment  at 
the  root  of  the  instrumented  blade  were  measured  by  a  five -component 
balance  system  which  held  the  root  of  the  blade.  Rotor  azimuthal  position, 
root  angle  of  attack  of  each  of  the  two  blades,  and  rotor  shaft  angle  (in 
the  plane  of  the  windstream)  were  remotely  controlled  from  a  console 
outside  the  tunnel. 

Figure  2  presents  a  dimensional  drawing  giving  the  pertinent 
dimensions  of  the  test  installation  in  the  wind  tunnel.  Figure  3  presents 
a  close-up  view  of  the  rotor-head  assembly  showing  the  various  drive 
systems.  The  blade  pitch  and  shaft  tilt  drives  were  linear  ball-screw 
actuators  that  were  designed  to  change  the  respective  angles  at  a  rate  of 
1/2  degree  per  second,  and  the  azimuth  drive  was  designed  to  rotate  the 
rotor  system  at  a  rate  of  1  degree  per  second.  The  blade  pitch-control 
actuators  were  driven  by  a  high-speed  motor,  and  thus  a  large-speed 
reduction  was  required  to  obtain  the  slow  rates  of  blade  pitch-angle 
change.  However,  because  of  the  extremely  small  amount  of  friction 
in  the  ballscrew  actuator,  the  use  of  a  brake  was  required  to  keep  the 
applied  aerodynamic  load  from  rotating  the  actuator  system  after  the 
drive  motor  had  been  turned  off.  This  was  accomplished  by  means 
of  a  small  friction  brake  that  was  controlled  through  a  solenoid  that 
released  the  brake  when  electrical  power  was  applied  to  the  drive 
motor.  As  can  be  seen  from  Figure  3,  the  blade  pitch  drive  for  the 
instrumented  blade  (designated  the  red  blade)  was  connected  to  the 
inner  housing  of  the  balance  system. 

Figure  4  shows  the  overall  details  of  the  balance  system.  The 
outer  balance  housing  was  fixed  to  the  head  assembly,  and  the  inner  bal¬ 
ance  housing  rotated  inside  the  outer  housing  on  two  large  radial  bearings. 
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One  side  of  each  balance  flexure  was  attached  to  the  inner  housing,  and 
the  other  side  of  the  flexure  was  attached  to  the  blade  grip  retention. 

The  center  of  rotation  of  the  balance  was  coincident  with  the  quarter 
chord  of  the  blade.  As  previously  noted,  the  blade  angle  of  attack  was 
changed  by  rotating  the  entire  blade-balance  unit  (inner  balance  housing). 
The  balance,  therefore,  always  measured  the  forces  and  moments  per¬ 
pendicular  and  parallel  to  the  root  chord  of  the  blade. 

The  disc  brake  shown  in  Figure  4  was  used  to  lock  out  the  small 
amount  of  blacklash  in  the  azimuth  drive  train  once  the  desired  rotor 
azimuth  position  had  been  obtained.  The  brake  was  an  electrohydraulic 
system  remotely  controlled  by  the  test  operator  in  the  control  room. 

Figure  5  is  a  dimensional  drawing  which  shows  the  arrangement 
of  the  force-beam  flexures  in  the  balance  system  and  a  detailed  sketch  of 
a  typical  force-beam  flexure.  As  can  be  seen  from  the  force-beam  lay¬ 
out,  the  root  edgewise  moments  and  the  root  beamwise  moments 

( B ma)  were  measured  by  pairs  of  force  beams  on  approximately  2-foot 
centers.  The  root  edgewise  force  (PPa)  was  measured  by  a  single  force 
beam,  and  the  root  normal  force  (nfp)  and  root  pitching  moment  (p/*a) 
were  both  measured  by  the  same  set  of  force  beams  spaced  on  16-inch 
centers.  To  obtain  the  APA,  the  strain  signals  from  the  two  force  beams 
were  added  electrically;  to  obtain  the  pitching  moment,  the  strain  signals 
were  subtracted  electrically. 

The  detailed  sketch  of  a  typical  force  beam  shows  how  a  built-in 
stress  concentration  was  used  to  facilitate  the  attainment  of  a  large  sig¬ 
nal  from  the  strain  gages  while  maintaining  a  very  rigid  balance  system. 
While  a  certain  amount  of  crosstalk  is  obtained  from  such  a  balance 
arrangement,  it  has  been  shown  from  experience  that,  by  means  of  a 
proper  calibration  procedure  whereby  all  cross  couplings  are  meas- 
sured,  very  accurate  force  and  moment  measurements  can  be  obtained. 
The  load  ranges  to  which  the  present  balance  system  was  designed  are 
as  follows: 


root  normal  force 
root  edgewise  force 
root  beamwise  moments 
root  edgewise  moments 
root  pitching  moments 


(*/pp) 

1  1,  200  lbs 

{&*) 

+ 500  lbs 

(an*) 

t  11, 000  ft  -lbs 

(*"*) 

+  5,  000  ft -lbs 

(p*p) 

±920  ft -lbs 

The  rotor  blades  that  were  used  for  the  tests  were  the  instru¬ 
mented  UH-1A  blades  previously  used  in  1962  to  obtain  the  blade  air¬ 
load  data  obtained  in  flight  and  presented  in  Reference  1.  A  complete 
and  detailed  description  of  the  rotor  blades  can  be  found  in  that  refer¬ 
ence.  The  blade  geometric  properties  are  presented  in  Figure  6,  the 
blade  mass  distribution  is  presented  in  Figure  7,  the  blade  mass  inertia 
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distribution  is  presented  in  Figure  8,  the  blade-bending  stiffness  dis¬ 
tributions  are  presented  in  Figure  9,  and  the  blade  torsional  stiffness 
distribution  is  presented  in  Figure  10.  In  the  configuration  with  which 
the  present  tests  were  conducted,  the  15-1 /4-inch  constant-chord  rotor 
blade  has  a  diameter  of  approximately  45  feet.  The  trim  tab  (similar 
to  the  one  shown  on  the  dummy  blade  in  Figure  1)  was  removed  from  the 
instrumented  blade  so  that,  when  the  blade's  trailing  edge  was  facing  the 
airstream,  the  flow  would  not  be  disturbed  by  the  protruding  tab. 

The  transducers  located  on  the  test  blade  were  of  two  types; 
strain  gages  to  measure  the  beamwise,  edgewise,  and  torsional  strains 
at  nine  spanwise  stations,  and  pressure  gages  to  measure  the  chordwise 
pressure  distributions  at  six  spanwise  stations.  Figure  11  shows  the 
spanwise  location  of  the  strain  gages  and  pressure  gages,  and  Table  I 
lists  the  chordwise  location  and  pressure  range  of  each  pressure  gage 
on  the  blade.  Figure  12  shows  both  the  pressure-gage  and  the  strain- 
gage  installation  at  approximately  85-percent  span.  As  can  be  seen, 
the  strain-gage  installation  is  a  relatively  standard  one  using  four  active 
gages  for  each  of  the  three  independent  strain  readings.  This  type  of 
strain-gage  installation  was  used  instead  of  rosettes,  as  the  majority 
of  the  strain-gage  instrumentation  previously  installed  for  use  in  the 
flight-loads  program  was  not  of  the  rosette  type  (Reference  1)  and  was 
adaptable  for  the  present  tests. 

In  the  pressure-gage  installations,  two  types  of  gages  were 
used:  a  NACA -developed  miniature  differential  pressure  transducer 
previously  installed  in  the  blade  (Reference  1)  and  an  absolute  pressure 
transducer  (Reference  2).  The  differential  gages  measured  the  differ¬ 
ences  between  the  upper  and  lower  surface  pressures.  The  absolute 
pressure  gages  were  used  to  measure  the  differences  between  the 
wind-on  and  wind-off  surface  pressures.  The  absolute  pressure 
transducers  were  installed  only  at  the  85-percent  span  station  in  an 
attempt  to  obtain  a  more  complete  definition  of  the  pressure  distri¬ 
bution  than  would  have  been  available  from  the  existing  differential 
gages.  As  can  be  seen  from  Figure  12,  a  chordwise  epoxy  fairing 
was  made  for  the  85 -percent  span  station  in  an  attempt  to  mount  the 
absolute  surface  gages  external  to  the  blade  without  disrupting  the 
airflow.  Prior  to  testing,  the  recesses  in  the  epoxy  fairing  for  the 
gages  were  filled  with  a  soft  plastic  to  smooth  the  contour. 

A  10-volt,  3,  000-cps  excitation  supplied  by  a  power  amplifier 
was  used  to  excite  both  the  strain  gages  and  pressure  transducers. 

Figure  13  shows  the  installation  of  the  transducers  used  to 
measure  the  blade  pitch  angles  of  the  instrumented  and  dummy  blades 
with  respect  to  the  rotor -head  baseplate.  The  transducers  for  both 
blades  were  one-turn,  precision  rotary  potentiometers  mounted  at 
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Table  I 

LOCATIONS  AND  RANGES  OF  PRESSURE  GAGES 


S PANWISE  STATION  (%  radius) 

CHORDWISE  STATION 

(%  chord; 

11.1 

(psi) 

55.7 

(psi) 

75.4 

(psi) 

85.2 

(psi) 

90.1 

(Psi) 

95.0 

(psi) 

2.1 

- 

15 

15 

15 

15 

15 

3.7 

8 

- 

- 

15 

- 

- 

9.0 

- 

15 

15 

15 

15 

15 

12.9 

- 

- 

- 

15 

- 

- 

16.9 

8 

8 

- 

15 

15 

15 

22.6 

8 

- 

15 

15 

- 

33.5 

4 

- 

- 

- 

- 

- 

117.8 

- 

- 

- 

4 

- 

- 

62.0 

1 

2 

2 

2 

4 

68  q  UPPER 
b8*3  SURFACE 

- 

- 

5-20* 

- 

- 

go  q  LOWER 
b8,3  SURFACE 

- 

- 

5-20* 

- 

- 

75.il 

- 

- 

• 

2 

_ 

82  0  UPPER 

8A0  SURFACE 

- 

- 

- 

5-20* 

- 

- 

85.9 

1 

- 

- 

- 

- 

- 

88.0 

- 

1 

1 

2 

2 

2 

Qii  1  UPPER 
SURFACE 

- 

- 

- 

5-20* 

1 

1 

qii  |  LOWER 

34,1  SURFACE 

- 

- 

- 

5-20* 

B 

■ 

97  0  UPPER 

3'-°  SURFACE 

- 

- 

- 

5-20* 

■ 

B 

97  0  LOWER 

3/-°  SURFACE 

- 

- 

- 

5-20* 

■ 

■ 
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the  center  of  rotation  of  the  inner  housing  of  the  balance  for  the 
instrumented  blade  and  at  the  center  of  rotation  of  the  blade-bearing 
housing  for  the  dummy  blade.  Both  potentiometers  were  powered  by 
a  stabilized  40 -volt  dc  power  unit. 

Figure  14  is  a  close-up  photograph  of  the  rotor  head  installed 
in  the  tunnel  prior  to  testing.  The  transducers  used  to  measure  the 
shaft  angle  a  id  azimuth  angle  can  be  clearly  seen  in  this  picture.  The 
transducer  used  to  measure  the  shaft  angle  was  a  270-degree  precision 
rotary  potentiometer  mounted  at  the  center  of  the  shaft  tilt  bearing  housing. 
The  azimuth  angle  transducer  was  a  10-turn  rotary  potentiometer  driven 
through  a  friction  contact  with  the  azimuth  brake  disc  (Figure  4).  One 
revolution  of  the  rotor  head  corresponded  to  approximately  eight  turns 
of  the  rotary  potentiometer. 

Since  the  azimuth  potentiometer  was  driven  by  the  brake  disc 
through  a  friction  contact,  a  wind  screen  (not  shown  in  the  photograph) 
was  used  to  prevent  the  airstream  from  altering  the  contact  pressure 
and  a  remote  means  (optical)  of  checking  on  possible  slippage  was  used 
to  monitor  the  performance  of  the  transducer.  A  fuller  discussion  of 
the  remote  optical  check  system  is  presented  in  the  section  describing 
the  operational  procedures  followed  during  the  wind  tunnel  tests. 

All  transducer  signals,  position  transducers,  pressure  gages, 
strain  gages,  etc.  ,  were  recorded  both  on  a  digital  readout  system 
installed  in  the  control  room  of  the  NASA -Ames  40  - by-80 -foot  full-scale 
wind  tunnel  (Datex  II)  and  by  two  50-channel,  type  7-119  oscillographs. 

The  two  readout  systems  were  employed  to  insure  that,  if  the  primary 
digital  data  recording  system  (Datex  II)  malfunctioned,  backup  direct 
analog  traces  would  be  available  to  obtain  the  desired  data.  In  addition, 
if  other  than  "static  conditions"  prevailed,  time  analog  traces  would  be 
recorded.  The  latter  was  the  primary  reason  that  the  analog  backup 
system  was  believed  to  be  eeded,  as  it  was  suspected  that  the  blade 
would  not  remain  stationary,  under  certain  test  conditions,  for  the 
relatively  long  time  required  by  the  Datex  II  unit  to  punch  all  the  data 
on  digital  cards  (approximately  12  seconds).  Because  of  mismatched 
impedances,  the  analog  and  digital  data  could  not  be  recorded  simul¬ 
taneously.  A  switching  system  was  developed,  therefore,  to  place  all 
the  transducer  signals  on  either  the  Datex  II  system  or  the  analog 
recording  system.  Figure  15  presents  a  block  diagram  of  the  entire 
instrumentation  system  and  indicates  the  switching  arrangement  that 
was  used. 

Figure  16,  a  photograph  of  the  instrumentation  installed  in  the 
wind  tunnel  control  room,  shows  the  eight  banks  of  amplifiers  and  the 
switching  box  for  signal  transfer  to  either  the  digital  recorder  or  the 
CEC  recording  oscillographs, which  are  also  shown  in  the  photograph. 

Figure  17  is  a  photograph  of  the  Datex  II  strip  chart  readout, 
which  permitted  an  immediate  inspection  of  the  data  that  were  being 
punched  for  each  channel  on  the  cards  during  every  cycle  of  data  taking. 
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Figure  18  is  a  photograph  of  the  test  operator's  control  panel. 
All  the  controls  and  appropriate  readouts  of  blade  pitch  angle  and  blade 
stress  for  the  dummy  blade  (white  blade)  were  on  the  left, and  the  control 
panel  and  readouts  for  the  instrumented  blade  (red  blade)  were  on  the 
right.  The  readouts  and  controls  for  azimuth  angle  and  shaft  tilt  angle 
common  to  both  blades  were  in  the  center  part  of  the  control  panel. 

Each  drive  system  had  two  switches  :  one  to  turn  the  power 
on  or  off  to  the  respective  control,  and  a  three-position, spring-loaded 
switch  to  drive  the  control  in  either  a  plus  or  a  minus  direction.  At  the 
top  center  of  the  control  panel  was  a  multichannel  selector  switch  which 
allowed  the  appropriate  signals  to  be  accurately  read  out  on  a  digital 
voltmeter  ,  The  blade  pitch  angles,  shaft  tilt  angle,  azimuth  angle, 
control  system  voltages,  and  instrumentation  voltages  were  all  set  and 
checked  by  means  of  the  digital  voltmeter  so  that  accurate  parameter 
values  could  be  obtained.  It  .  hould  be  noted  that  the  readout  dials  on 
the  control  panel  were  used  only  as  a  quick  visual  indication  of  the  rotor 
and  blade  positions  and  that  the  digital  voltmeter  and  the  recording 
systems  were  used  to  set  and  record  accurate  values  of  the  parameters. 

B.  EQUIPMENT  CHECKOUT  AND  CALIBRATION 


Prior  to  shipment  of  the  rotor  system  to  the  NASA  -Ames  40- 
by-80-foot  wind  tunnel  for  tests,  the  entire  system  was  assembled  and 
checked  out  for  proper  operation.  In  addition,  all  major  instrumen¬ 
tation  systems  were  completely  and  carefully  calibrated  so  that  only 
check  calibration  would  be  required  after  installation  in  the  tunnel. 

The  primary  measurement  subsystems  were  calibrated  independently, 
and  then  as  a  complete  unit.  These  subsystems  were  the  amplifier - 
analog  recording  system,  the  force-moment  balance  at  the  root  of  the 
instrumented  blade,  the  various  position  transducers,  the  pressure 
transducers  on  the  blade,  and  the  strain  gages  mounted  on  the  instru¬ 
mented  blade.  The  digital  recording  system  installed  in  the  control 
room  of  the  wind  tunnel  at  Ames  was  set  and  calibrated  after  the  rotor 
system  had  been  mounted  in  the  test  section.  The  following  sections 
of  the  report  will  discuss  the  calibration  procedures  and  results  obtained 
for  each  of  the  measurement  subsystems  that  were  calibrated  independently. 

Amplifier -Analog  Recording  System 

The  amplifier -recording  system  was  calibrated  in  two  steps. 

The  amplifier  input -to -output  voltage  amplification  for  each  attenuator 
setting  was  checked  to  insure  that  the  specified  ratios  were  accurate  to 
within  1  percent, and  then  the  galvometer  deflections  in  the  recording 
oscillograph  were  determined  for  each  transducer  channel  in  terms  of 
the  voltage  output  of  the  amplifiers.  Table  II  presents  the  recorded 
trace  deflection  of  each  transducer  channel  per  volt  input  to  the  asso¬ 
ciated  amplifier  set  at  an  attenuation  of  unity. 
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Table  II 

AMPLIFIER-OSCILLOGRAPH  CALIBRATION  CONSTANTS 


Table  II  -  Continued 


The  Balance  System 


The  balance  system  was  calibrated  by  the  instrumentation- 
calibration  group  of  CAL's  Transonic  Wind  Tunnel  Department,  which 
has  been  responsible  for  the  instrumentation  and  calibration  of  all  the 
balance  systems  developed  for  use  in  that  facility.  The  balance  was 
calibrated  by  applying  known  loads  in  each  of  the  six  balance-component 
directions  independently.  The  balance  was  also  subjected  to  known  com¬ 
bined  loadings.  The  data  from  these  calibrations  were  used  to  determine 
the  matrix  relating  the  applied  loads  to  the  balance  output  signals.  The 
differences  in  the  calibration  constants  between  positive  and  negative 
loadings  were  also  determined.  The  balance  conversion  matrix  (balance 
beam  transducer  outputs  to  loads)  as  determined  during  calibration  is 
presented  in  Figure  19.  The  plus  or  minus  sign  in  front  of  each  row 
of  the  conversion  matrix  indicates  which  set  of  numbers  in  each  row 
should  be  used  for  either  a  positive  or  a  negative  load.  Care  had  to  be 
used  during  data  reduction,  therefore,  to  insure  that  the  correct  set 
of  matrix  elements  was  used.  It  is  noted,  however,  that  if  the 
balance  beam  signals  Xf  through  Xt  resulting  from  a  positive  normal 
force  and  drag  force  were  used  with  the  matrix  elements  for  negative 
forces,  positive  forces  would  be  obtained  although  their  magnitudes 
would  be  incorrect.  The  results  would  indicate,  however,  that  the 
matrix  elements  for  a  positive  force  should  have  been  used  to  obtain 
the  correct  forces  and  moments.  The  automatic  digital  data  reduction 
program  always  used  the  coefficients  for  positive  forces  first  and  then 
the  sign  of  the  calculated  forces  was  checked  with  the  assumption.  If 
either  of  the  forces  was  found  to  be  negative,  the  program  would  auto¬ 
matically  recalculate  the  loads  based  upon  the  other  set  of  coefficients. 

It  is  noted  from  inspection  of  the  balance  conversion  matrix 
(Figure  19)  that  it  is  not  diagonal  and  that  reasonably  large  interaction 
conversion  constants  are  present. 

At  the  completion  of  the  balance  calibration,  values  of  dummy 
load  resistors  that  could  be  applied  across  the  arms  of  each  balance 
beam  bridge  were  determined  for  full-scale  loads  so  that  electrical 
calibration  checks  of  the  balance  system  could  be  easily  made  during 
the  tests.  The  dummy  load  resistors  (DLR's)  that  were  established 
for  the  100-percent  load  condition  of  each  balance  beam  are  as  follows; 

Xt  88,  000  ohms  xf  36,  000  ohms 

X2  78,  000  ohms  xf  64,  000  ohms 

Xj  36,  000  ohms  40,  000  ohms 

The  calibration  of  the  balance  system  as  an  independent  unit 
indicated  repeatedly  that,  through  the  use  of  the  balance  conversion 
matrix,  the  forces  and  moments  could  be  determined  within  1  percent 
of  their  designed  full-scale  values.  Checks  made  during  calibration 
of  the  entire  system  established  that  the  loads  could  be  measured  to 
within  only  1-1/2  percent  of  the  designed  full-scale  values.  Using  this 
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latter  percentage  as  the  possible  error  range,  the  various  forces  and 
moments  have  the  following  error  bands: 


root  normal  force 

(NF*.) 

+  20  lbs 

root  edgewise  force 

(£FR) 

+  7  lbs 

root  beamwise  moments 

(Bhr) 

tl65  ft-lbs 

root  edgewise  moments 

(eMR) 

+75  ft-lbs 

root  pitching  moments 

( PMR.) 

+  15  ft-lbs 

Position  Transducers  Calibration 


The  azimuth  angle  transducer  was  calibrated  by  setting  the 
blade  at  prespecified  azimuthal  positions  and  recording  the  voltage 
output  from  the  transducer.  The  prespecified  azimuthal  positions 
consisted  of  twenty-four  15-degree  increments.  The  15-degree  incre¬ 
ments  were  marked  on  the  circular  rotor-head  baseplate.  This  base¬ 
plate  rotated  with  the  blade.  A  fixed  indicator,  having  a  scale  scribed 
on  it,  was  mounted  so  that  the  indicator  scale  and  baseplate  scale  were 
in  close  proximity  and  could  be  read  simultaneously.  The  indicator 
scale  was  divided  into  increments  of  0.  1  degree.  Thus,  the  prespeci¬ 
fied  azimuthal  positions  could  be  read  to  within  t  0,  1  degree.  The 
results  of  the  calibration  are  presented  in  Figure  20.  It  is  believed 
that,  on  the  basis  of  the  results  obtained  during  the  calibration,  the 
azimuth  angle  could  be  set  to  within  +  0.  5  degree  when  the  transducer 
output  was  read  on  the  digital  voltmeter  located  on  the  operator's 
console. 


The  transducer  for  the  shaft  tilt  angle,  as,  the  instrumented 
(red)  blade  pitch  angle,  and  the  dummy  (white)  blade  pitch  angle  were 
calibrated  using  an  inclinometer  capable  of  measuring  to  0.01  degree. 
The  results  of  these  calibrations  are  shown  in  Figures  21  through  23. 

It  is  believed,  on  the  basis  of  the  calibration,  that  the  shaft  tilt  angle 
and  blade  pitch  angles  were  measured  to  within  t  0.  20  degree  if  the 
linearized  curves  presented  in  Figures  21  through  23  are  used.  It 
is  noted  that  the  accuracy  of  determining  each  of  the  respective  angles 
was  significantly  improved  (+  0.  10  degree)  when  the  digital  voltmeter 
was  used  to  set  given  angles  for  which  a  specific  output  voltage  was 
known. 

Pressure  Gage  Calibration 

The  equipment  used  during  calibration  of  the  pressure  gages 
is  shown  in  Figure  24.  A  pressure  fitting  was  contoured  to  the  cross 
section  of  the  blade  so  that  one  whole  chordwise  station  could  be  cali¬ 
brated  at  one  time.  Since  a  vacuum  was  applied  to  the  gages,  a  good 
seal  between  the  pressure  fitting  and  the  blade  could  be  obtained  through 
the  use  of  a  fibrous  putty.  The  vacuum  gage  was  capable  of  indicating 
0.01  psi.and  0.  005  psi  could  be  easily  estimated.  All  gages  were 
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calibrated  at  0.  20-psi  increments  for  increasing  and  decreasing  pres¬ 
sures  up  to  a  maximum  pressure  of  4  psi.  The  results  of  those  cali¬ 
brations  for  typical  1  -,  2-,  4-,  8-,  and  15-psi  differential  gages  are 
presented  in  Figures  25  through  29;  for  the  5-  to  20-psi  absolute  gages, 
in  Figure  30.  Table  III  presents  the  calibration  constants  in  terms  of 
volts  per  psi  for  all  the  gages  mounted  on  the  blade. 

It  is  noted  from  the  results  presented  in  Figures  25  through 
30  that  the  transducer  outputs  were  surprisingly  linear  throughout  the 
pressure  range  of  interest.  It  is  noted,  for  example,  that  the  1-psi 
differential  gage  was  linear  up  to  three  times  its  designed  output,  and 
the  15-psi  gage  was  surprisingly  accurate  in  a  pressure  range  of  fl  psi. 
Although  it  has  not  been  indicated  in  the  figures,  the  repeatability  of 
each  of  the  gages  was  remarkably  good,  with  the  differences  being  too 
small  to  be  plotted.  On  the  basis  of  the  results  obtained,  it  is  believed 
that  the  pressure  gages  were  sufficiently  well  calibrated  so  that  the 
readings  for  all  gages  can  be  expected  to  have  an  error  of  no  more 
than  +0.  2  psi  for  the  pressure  range  of  interest  (0  to  +2.  5  psi). 

Calibration  of  Blade-Mounted  Strain  Gages 


The  calibration  constants  for  the  strain  gage  sets  which  meas- 
sured  the  beamwise,  edgewise,  and  torsional  moments  are  presented  in 
Figures  31  through  33.  These  constants  were  obtained  as  follows:  con¬ 
sider  the  beamwise  gage  set  calibration  constants  given  in  Figure  31  as 
typical.  The  beam  was  loaded  at  a  given  station  with  a  beamwise  moment. 
The  signal  from  the  gage  set  which  was  designated  to  measure  beamwise 
moments  was  recorded  for  each  instrumented  spanwise  station.  The 
signal  from  each  beamwise  gage  set  was  recorded  as  the  moment  was 
increased  at  the  given  station.  A  plot  of  signal  output  versus  moment 
applied  was  made  for  each  spanwise  station,  and  the  resultant  slope 
(volts /ft-lb)  was  determined.  The  point  of  application  of  the  moment 
was  changed  and  the  process  repeated.  The  slopes  obtained  at  each 
spanwise  station  for  the  various  load  stations  are  presented  in  Figure  31. 
For  purposes  of  data  reduction,  the  slopes  at  each  station  were  averaged; 
the  averaged  values  are  presented  in  Table  IV.  Coupling  between  the 
various  moments  was  neglected. 

To  determine  the  moments  induced  by  the  blade  weight,  the 
blade  was  supported  in  a  horizontal  position.  This  position  was  used 
as  the  zero  moment  reference  position  during  pretest  calibrations. 

The  moments  due  to  the  blade  weight  were  measured  with  respect  to 
this  reference  position.  The  beamwise  and  edgewise  deflections  of 
the  beam  were  also  measured  and  are  presented  in  Figures  34  and  35, 
respectively.  The  torsional  deflection  of  the  blade  due  to  its  weight 
was  not  measurable. 

During  the  tunnel  tests,  however,  the  zero  reference  position, 
for  reasons  of  convenience,  was  chosen  as  that  of  the  blade  deformed 
under  its  own  weight.  The  moment  data  were  then  appropriately  modified 
to  account  for  the  weight  contribution. 
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Table  UL 

PRESSURE- GAGE  CALIBRATION  CONSTANTS 


L 

($  SPAN) 

GAGE 

OCATION 

CHORD) 

CALIBRATION 

CONSTANT 

,  VOLTS/PS  1  , 

(ATTENUATIONS) 

85.2 

47.8 

0.275 

85.2 

62.0 

0.620 

85.2 

68.9 

UPPER 

SURFACE 

0.060 

85.2 

68.9 

LOWER 

SURFACE 

0.069 

85.2 

75.4 

0.520 

85.2 

82.0 

UPPER 

SURFACE 

0.051 

85.2 

88.0 

0.326 

85.2 

94.1 

UPPER 

SURFACE 

0.055 

85.2 

94.1 

LOWER 

SURFACE 

0.076 

85.2 

97.0 

UPPER 

SURFACE 

0.044 

85.2 

97.0 

LOWER 

SURFACE 

0.053 

90.1 

2.1 

0.163 

90.1 

9.0 

0.123 

90.1 

16.9 

0.113 

90.1 

22.6 

0.165 

90.1 

62.0 

0.334 

90.1 

88.0 

0.435 

95.0 

2.1 

0.154 

95.0 

9.0 

0.160 

95.0 

16.9 

0.075 

95.0 

88.0 

0.530 

LOCATION 


{%  SPAN)  CHORD) 


CALIBRATION 

CONSTANT 

VOLTS/PS  I 
ATTEHUAT I 0N=1 ) 


41.1 

3.7 

41.1 

16.9 

41.1 

33.5 

41.1 

62.0 

41.1 

85.9 

55.7 

2.1 

55.7 

3.7 

55.7 

16.9 

55.7 

22.6 

55.7 

62.0 

55.7 

88.0 

75.4 

2.1 

75.4 

9.0 

75.4 

62.0 

75.4 

88.0 

85.2 

2.1 

85.2 

3.7 

85.2 

9.0 

85.2 

12.9 

85.2 

16.9 

85.2 

22.6 

Table  II 

STRAIN-GAGE  CONVERSION  CONSTANTS 


$  SPAN 

BEAMWISE 

(volts/ft -lb  X  I03) 

EDGEWISE 

(vol t* /ft- lb  x  I03) 

TORSION 
(vol  tt/f  t  -lb 

93.0 

1.560 

0.097 

1.296 

8b. 9 

1.790 

0.182 

1.366 

80.7 

1.340 

0.239 

1.332 

70.8 

1.819 

0.267 

1.296 

61.0 

1.670 

0.280 

1  .M2 

51.1 

1.760 

0.299 

1.230 

37.3 

1.890 

0.200 

1 .  MO 

28.il 

1.810 

0.197 

1.262 

16.6 

0.6*10 

0. 1 49 

0.197 

It  should  be  noted  at  this  point  that  a  considerable  effort  was 
expended  in  an  attempt  to  obtain  blade  load  distribution  data  from  a 
knowledge  of  the  blade  deflection  data  as  determined  from  the  strain 
gages.  This  effort  failed  to  yield  any  meaningful  data.  It  is  believed 
that  the  reasons  for  the  failure  were  twofold. 

The  first  reason  was  the  mathematical  nature  of  the  equations 
describing  the  relationship  between  the  loads  and  moments  and  deflec¬ 
tions.  This  relationship  can  be  written  in  matrix  form  as 

l‘l  =  H  W 

where 

{£]  column  matrix  of  blade  deflections 

[l]  column  matrix  o:  blade  loads  and  moments 

[e]  square  matrix  of  blade  influence  coefficients 


The  matrix  [c]  can  be  determined  experimentally.  Mathematically,  it 
is  possible  to  obtain  the  inverse  of  the  matrix[c]  and,  hence,  to  obtain 
load  information  from  deflection  data  (i.  e. ,  strain  signals).  However, 
practical  difficulties  arise  in  obtaining  the  inverse  of  the  matrix  [c]  when 
[c]  has  off-diagonal  elements  which  are  large  compared  to  the  main  diag¬ 
onal  elements.  The  matrix  fc]  is  then  said  to  be  ill-condi.tioned,and  an 
accurate  inverse  is  obtainable  only  under  very  stringent  conditions.  For 
the  case  in  hand,  these  conditions  were  not  met.  A  more  detailed  dis¬ 
cussion  of  this  problem  is  presented  in  Reference  18. 

A  second  reason  for  the  lack  of  success  of  the  effort  was  that, 
for  the  stopped-rotor  configurations,  blade  deflections  could  arise 
which  result  in  nonlinear  relations  between  the  applied  loads.  Since 
the  entire  approach  is  based  upon  the  theory  of  linear  superposition, 
it  did  not  succeed.  As  an  example,  consider  the  pitching  deflection 
of  the  blade  root.  A  contribution  to  the  pitching  deflection  from  the 
drag  load  arises  when  there  is  a  spanwise  bending  deflection.  However, 
this  spanwise  deflection  is  due  principally  to  the  blade  normal  force. 
Hence,  a  pitching  deflection  arises  which  is  proportional  to  the  product 
of  the  drag  and  normal  force  acting  on  the  blade. 

C.  WIND  TUNNEL  TEST  PROCEDURES 


After  the  test  rotor  had  been  mounted  in  the  wind  tunnel,  final 
calibrations  of  the  shaft  tilt,  blade  pitch,  and  azimuth  angle  transducers 
were  conducted.  Check  calibrations  were  also  conducted  on  all  trans¬ 
ducers  to  insure  that  they  were  operating  correctly  and  that  their  cali¬ 
bration  constants  had  not  changed. 
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Since  the  axis  of  the  force  balance  at  the  root  of  the  instru¬ 
mented  blade  rotated  with  the  blade  during  pitch  changes,  the  gravity 
effect  of  the  blade  on  the  balance  had  to  be  determined  over  the  range 
of  blade  pitch  angles,  shaft  tilt  angles,  and  azimuth  angles  at  which 
the  tests  were  to  be  conducted.  A  series  of  "static  runs"  (no  free- 
stream  velocity)  were  conducted  prior  to  the  "wind-on"  tests  to  obtain 
the  balance  tares  for  the  expected  range  of  rotor  parameters.  Over 
a  hundred  data  points  were  taken  during  the  "static  runs"  to  provide 
balance  tare  values  for  various  combinations  of  ihe  rotor  control 
angles,  as  ,  &K  ,  * . 

Prior  to  each  data  run  (a  data  run  is  measured  from  the 
time  the  air  is  turned  on  to  the  time  the  tunnel  is  shut  down),  the 
instrumented  blade  was  put  at  zero  azimuth  angle, and  the  shaft  tilt 
angle  and  blade  root  angle  were  also  placed  at  zero.  All  instrumen¬ 
tation  channels  were  balanced, and  all  transducer  signals  were  then 
recorded  for  a  standard  attenuation  on  each  channel.  These  signals 
were  considered  to  be  initial  zeros.  At  the  end  of  each  data  run,  the 
same  procedure  was  carried  out  to  obtain  final  zeros,  except  the  instru¬ 
mentation  was  not  rebalanced.  Comparison  of  before-and-after  run 
zeros  permitted  "instrumentation  drift"  to  be  evaluated.  It  was  found 
to  be  negligible  in  all  cases. 

The  tests  were  conducted  by  an  operator  who  sat  at  the  control 
console  in  the  forward  portion  of  the  run  shack,  where  he  could  have  a 
full  view  of  the  rotor.  Figure  36  is  a  photograph  of  the  test  operator's 
station.  The  transit  to  the  right  of  the  test  operator's  station  was  used 
to  check  for  slippage  of  the  friction  wheel  on  the  azimuth  transducer 
during  a  run.  This  was  accomplished  by  setting  the  rotor  at  ^  =  90 
degrees  with  the  digital  readout  and  then  viewing  two  scribe  marks  on 
the  rotor  head  (one  on  the  rotating  side  and  one  on  the  nonrotating  side) 
to  see  if  they  were  properly  aligned.  If  the  offset  between  the  scribe 
marks  indicated  that  more  than  approximately  1/4  degree  of  slippage 
had  occurred,  the  tunnel  was  shut  down  and  the  transducer  reset.  This 
had  to  be  done  only  once  during  the  tests,  when  the  amount  of  slippage 
indicated  that  the  azimuth  angle  was  approximately  0.  4  degree  in  error. 

The  test  operator  communicated  with  the  instrumentation  and 
data  engineers  located  in  the  back  of  the  run  shack  (Figure  16)  through 
a  sound -powered  telephone  set 

The  test  procedure  used  to  obtain  the  data  varied  because  of 
limits  imposed  by  static  loads  and  instabilities.  To  obtain  the  data  at 
the  low  dynamic  pressure,  ^  =  15  psf,  the  red  blade  pitch  angle  and  the 
shaft  tilt  angle  were  held  fixed  and  the  azimuth  angle  was  varied.  For 
nominal  dynamic  pressures  of  30  and  45  psf,  the  azimuth  angle  was 
held  fixed  and  the  red  blade  pitch  angle  and  shaft  tilt  angle  were  varied 
to  obtain  the  data  at  the  desired  conditions. 
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The  primary  limits  that  were  imposed  on  the  conditions  at 
which  data  could  be  obtained  were  dynamic  instabilities  encountered 
in  the  azimuth -angle  range  of  225  to  300  degrees,  and  a  prespecified 
maximum  beamwise  stress  of  30,  000  psi  at  the  blade  critical  section. 

The  blade  critical  station  was  chosen  to  be  the  28.  4-percent  span  since 
this  was  the  inboard-most  station  prior  to  the  increase  of  blade  stiffness 
at  the  root.  Thus,  the  largest  blade  stresses  would  probably  occur  here. 
The  maximum  spanwise  stress  was  arrived  at  by  taking  three-quarters 
of  the  allowable  stress  for  the  weakest  material  used  in  the  blade.  The 
moment  at  the  critical  station  which  would  cause  a  stress  of  30,  000  psi 
was  approximately  5,  500  ft-lbs. 

The  nominal  conditions  for  which  data  points  were  obtained  at 
nominal  dynamic  pressures  of  15,  30  and  45  psf  are  presented  in  Tables 
V,  VI,  and  VII,  respectively.  In  all,  over  700  data  points  were  obtained 
during  the  tests. 

D.  PRESENTATION  AND  DISCUSSION  OF  RESULTS 


The  results  that  will  be  presented  and  discussed  in  this  report 
are  primarily  those  obtained  from  the  balance  measurements.  All  the 
basic  data  that  were  obtained  during  the  test  program  will  be  held  on 
file  at  USAAVLABS.  Figure  37  presents  a  typical  printout  of  the  various 
data  items  that  were  obtained  at  each  test  condition.  As  noted  in  Fig¬ 
ure  37,  the  test  parameters  are  listed  in  the  upper  right-hand  corner 
of  the  data  page.  These  parameters  consist  of  the  dynamic  pressure 
(lbs/ft^),  blade  azimuth  position  (degrees),  shaft  tilt  angle  (degrees), 
blade  root  pitch  angle  (degrees),  and  the  run  and  point  numbers  which 
define  the  sequence  of  the  actual  test  points.  Next  is  presented  the 
matrix  of  pressures,  differential  and  absolute,  in  psi.  Each  column 
of  the  matrix  represents  the  chordwise  pressure  distribution  at  a  given 
spanwise  station.  Most  of  the  pressures  presented  are  differential; 
however,  those  denoted  by  the  symbols  TS  and  LS  are  absolute  pres¬ 
sures.  The  TS  symbol  indicates  the  pressure  measured  on  the  upper 
surface  of  the  blade  at  the  given  station,  and  the  LS  denotes  the  lower 
surface  pressure.  Positive  readings  denote  pressures  which  correspond 
to  positive  lift;  therefore,  pressures  below  atmospheric  are  positive  for 
the  upper  surface  while  pressures  above  atmospheric  are  positive  on  the 
lower  surface. 

Considerable  care  should  be  exercised  in  interpreting  the 
pressure  distributions  that  were  obtained,  particularly  for  azimuth 
positions  ranging  from  180  degrees  to  330  degrees.  This  caution  is 
given  as  the  trailing  edge  of  the  blade  was  blunted  because  of  the  signif¬ 
icant  number  of  gage  instrumentation  wires  that  were  bonded  along  the 
span  near  the  trailing  edge.  The  effect  of  the  blunted  trailing  edge 
might  be  significant  when  the  flow  is  such  as  to  make  the  trailing  edge 
a  leading  edge.  Furthermore,  upon  inspection  of  the  data,  it  was  found 
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Table  1 

TABULATION  OF  NOMINAL  TEST  CONDITIONS 
(dynamic  PRESSURE,  J  -  15  psf) 

X... INDICATES  A  TEST  POINT 

AZIMUTHAL 

POSITION 

(/  ) 
deg 

SHAFT 

ANGLE 

(fs) 

deg 

BLADE  ROOT 

ANGLE,  & 

,  deg 

-3 

D 
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+7 
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X 
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0 

X 

X 

X 

X 
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-5 

-10 

n 

120 
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-10 

H 

150 

+10 

■ 

X 

+5 

X 

0 

I 

X 

X 

X 

X 

X 

-5 
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n 
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-10 

H 
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+10 
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+5 
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X 

X 

X 
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X 

-10 

D 
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+10 

■ 
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X 

X 

X 

X 

X 

-5 

-10 

n 

270 

+10 

■ 

X 

+5 

X 

0 

i 

X 

X 
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X 

X 

-5 

-10 

a 
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Table  JILL' 

TABULATION  OF  NOMINAL  TEST  CONDITIONS 

(DYNAMIC  PRESSURE,  <f.  -  46  psf) 

X... INDICATES  A  TEST  POINT 


that  very  low  pressure  levels  were  indicated  in  certain  azimuthal  posi¬ 
tions.  Consequently,  great  reliance  was  placed  on  the  integrated  loads 
measured  by  the  balance  system.  Analyses  of  the  pressure  data  were 
utilized  primarily  to  assist  in  the  interpretation  of  the  other  measure¬ 
ments  wherever  possible. 

The  spanwisc  distributions  of  the  blade  flatwise,  edgewise, 
and  torsional  structural  moments  are  presented  immediately  below  the 
pressure  matrix, as  shown  in  Figure  37.  As  previously  noted,  these 
moments  do  not  include  the  gravity  moments  due  to  blade  weight  ; 
thus,  they  are  the  true  aerodynamic  moments.  Finally,  the  balance 
data  are  presented.  The  balance  was  fixed  to  the  blade  root,  and  the 
data  presented  here  are  in  the  balance  frame  of  reference. 

It  was  necessary  to  correct  the  balance  readings  for  the  tares 
due  to  blade  weight.  Consideration  of  these  tares  resulted  in  a  modi¬ 
fication  of  the  balance  error  limits  from  those  presented  on  page  11 
(note  particularly  the  root  edgewise  force  and  root  edgewise  moment). 
The  balance  data  presented  in  the  tabular  listings  and  in  the  plots  to 
be  presented  herein,  therefore,  have  the  following  error  lirnits: 


root  normal  force 
root  edgewise  force 
root  beamwise  moments 
root  edgewise  moments 
root  pitching  moments 


(AW) 

+  25  lbs 

(bfr) 

+  15  lbs 

(8**) 

1 175  ft  -lbs 

(*"*) 

+  125  ft -lbs 

(PHR) 

+  1  5  ft  -  lbs 

Data  obtained  at  azimuth  angles  of  90,  270, and  180  degrees 
will  be  discussed  first.  For  azimuth  angles  of  90  and  266  degrees,  the 
various  forces  and  moments  as  measured  by  the  balance  system  have 
been  combined  to  yield  the  lift,  drag,  rolling  moment,  drag  moment, 
and  pitching  moment.  For  the  azimuth  angle  of  180  degrees,  only  the 
balance  measurements  are  presented. 

Figure  38  presents  the  variation  of  lift  with  blade  root  geo¬ 
metric  angle  of  attack  for  the  blade  at  azimuth  angles  of  90  and  266 
degrees.  The  sign  convention  of  the  composite  angle  of  attack  which 
is  made  up  of  the  blade  root  pitch  angle,  6*  ,  and  the  shaft  angle,  cr,  , 
for  both  azimuth  positions  is  indicated  on  the  respective  plots.  It 
should  be  remembered  that  this  blade  had  11  degrees  of  linear  twist 
(washout)  from  root  to  tip.  It  is  noted  that  the  data  presented  on  these 
plots  were  not  obtained  from  a  single  run  but  were  obtained  from  runs 
conducted  throughout  the  entire  series  of  wind  tunnel  tests  for  various 
combinations  of  shaft  angle  and  blade  pitch  angle.  It  was  gratifying 
to  observe  that  all  the  data  seemed  to  plot  on  fairly  smooth  curves 
with  only  the  anticipated  amount  of  scatter  present. 
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In  analyzing  the  results  as  presented  in  the  data  plotted  in 
Figure  38,  it  is  noted  that  the  apparent  integrated  nonlinear  stall 
effects,  as  measured  by  the  balance,  are  much  more  pronounced  at 
90  degrees  than  at  266  degrees,  in  that,  at  90  degrees,  the  lift  drop-off 
is  more  severe.  It  is  interesting  also  to  note  that  apparent  stall  occurs 
at  a  lower  root  geometric  angle  of  attack  when  the  trailing  edge  is  into 
the  airstream  than  when  the  normal  leading  edge  is  in  this  orientation. 
This  result  is  what  might  be  expected  based  upon  sectional  data  pre¬ 
sented  in  Reference  8.  The  nonlinear  stall  effects,  as  indicated  by 
the  data  presented  in  Figure  38  for  both  azimuth  positions,  are  more 
like  those  of  a  "soft  stall",  in  that  a  sharp  break  in  the  lift  curve  is 
not  obtained.  It  is  believed  that  this  "soft  stall"  characteristic  is  due 
to  the  fact  that  the  lift  is  the  air  load  integrated  over  the  whole  span; 
thus,  the  blade  is  effectively  slowly  stalling  with  increasing  angle  of 
attack. 


It  is  noted  that,  at  0  degrees  geometric  angle  of  attack,  the 
blade  is  developing  approximately  twice  the  lift  magnitude  at  an  azimuth 
position  of  266  degrees  than  it  is  at  90  degrees.  It  is  believed  that  this 
difference  is  due,  in  part,  to  the  blade  twist  resulting  from  the  aero¬ 
dynamic  pitching  moment.  At  90  degrees  azimuth,  the  aerodynamic 
pitching  moment  about  the  quarter  chord  at  0  degrees  geometric  angle 
of  attack  is  very  small;  while  at  266  degrees,  the  aerodynamic  pitching 
moment  is  significant  and  in  a  direction  so  as  to  increase  the  effective 
blade  angle  of  attack.  Calculations  based  on  the  measured  torsional 
strains  indicated  that  the  blade  tip  was  twisted  approximately  2  degrees, 
which  would  increase  the  total  blade  lift  at  y  =  266  degrees  by  approxi¬ 
mately  100  pounds.  This  increment  in  lift  due  to  twist  reduces  the 
differential  lift  magnitude  at  0  degrees  angle  of  attack  to  150  pounds 
between  the  90-  and  266-degree  azimuth  positions. 

It  was  suspected  that  this  remaining  difference  in  lift  might  be 
due  to  an  unknown  tunnel  flow  misalignment.  This  suspicion  was  sup¬ 
ported  by  the  fact  that  the  root  geometric  angles  of  zero  lift  differ  by 
3  degrees  for  the  90-  and  266-degree  azimuth  positions.  If  the  blade 
did  not  have  any  tip  effects  and  was  rigid,  the  root  angle  for  zero  lift 
for  both  azimuth  positions  would  be  one-half  of  the  blade  twist  or 
approximately  5-1/2  degrees.  This  angle  of  zero  lift  is  1-1/2  degrees 
higher  for  the  90 -degree  orientation  and  1-1/2  degrees  less  than  that 
measured  at  an  azimuth  angle  of  266  degrees.  While  flow  angularity 
could  explain  the  noted  differences,  a  flow  angularity  of  this  magnitude 
could  not  be  justified  on  the  basis  of  (low  calibrations  conducted  pre¬ 
viously  by  NASA -Ames. 

The  total  aerodynamic  rolling  moments  measured  by  the 
balance  at  ?  -  90  and  266  degrees  are  shown  in  Figure  39.  The  inter¬ 
esting  point  to  be  made  concerning  these  data  is  that  the  angle  for  zero 
rolling  moment  is  almost  the  same  as  that  for  zero  lift  at  f  =  266  degrees, 
and  only  1  degree  different  from  that  for  zero  lift  at  f  =  90  degrees.  The 
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difference  in  the  angle  of  zero  lift  and  zero  rolling  moment  would  be 
expected  to  be  approximately  2  degrees  for  a  rigid  blade  having  a  linear 
twist  distribution  of  1 1  degrees.  It  is  suspected  that  the  disagreements 
might  be  attributed  to  three-dimensional  elastic  effects. 

The  pitching  moments  about  the  quarter  chord  versus  root 
geometric  angle  of  attack  are  presented  for  azimuth  angles  of  90  and 
266  degrees  in  Figure  40.  It  is  noted  that  the  pitching  moments  at 
=  266  degrees  are  much  larger  than  at  f  -  90  degrees.  This  is  as 
expected  because  of  the  significant  difference  in  the  distance  between 
the  center  of  pressure  and  the  quarter  chord  at  f  =  90  and  266  degrees. 
The  difference  in  the  variation  of  the  pitching  moment  with  root  geo¬ 
metric  angle  of  attack  is  also  significant  between  f  =  90  and  266  degrees. 
Sectional  data  presented  in  Reference  8,  however,  tend  to  confirm  the 
differences  in  the  characteristics  that  were  measured.  The  large 
decrease  in  the  pitching  moment  above  13  degrees  at  ^  =90  degrees 
is  possibly  due  to  stall, which  tends  to  move  the  integrated  effective 
center  of  pressure  aft  of  the  quarter  chord  as  stall  progresses  along 
the  span  of  the  blade.  At  y/  =  266  degrees,  the  variation  of  the  pitching 
moment  with  angle  of  attack  would  be  expected  to  be  more  linear.  In 
this  case,  the  chordwise  shift  of  the  center  of  pressure  with  angle  of 
attack  is  a  smaller  percentage  of  the  distance  between  the  center  of 
pressure  and  the  quarter  chord  than  it  is  at  ^  =  90  degrees.  Thus,  a 
significant  nonlinearity  in  the  pitching  moment  about  the  quarter  chord 
would  not  occur  until  a  much  larger  portion  of  the  blade  has  become 
stalled.  In  addition,  it  is  noted  that  the  pitching  moment  at  ^  =  266 
degrees  would  not  change  sign  as  it  does  at  f  =  90  degrees,  as  the  blade 
stalls  with  increasing  angle  of  attack  since  the  center  of  pressure  will 
move  forward  only  to  approximately  the  50-percent  chord  station. 

The  variations  of  the  drag  force  with  blade  root  angle  are 
shown  in  Figure  41  for  V  =  90  and  266  degrees.  It  is  interesting  to 
note  that  the  root  geometric  angle  of  attack  for  minimum  drag  is 
approximately  the  same  for  both  azimuth  positions  and  that  both 
curves  are  approximately  symmetrical  about  the  angle  of  minimum 
drag.  It  is  somewhat  surprising,  based  on  sectional  data  presented 
in  Reference  8,  that  the  drag  force  at  f  =  266  degrees  was  not  signifi¬ 
cantly  larger  than  that  at  90  degrees.  On  the  basis  of  sectional  data, 
it  would  be  expected  that  the  drag  force  would  be  at  least  twice  as  large 
at  -  266  degrees  as  it  is  at  f  -  90  degrees  over  the  range  of  angle  of 
attack  for  which  results  were  obtained.  A  rational  explanation  for  the 
approximate  equivalence  of  the  drag  forces  at  the  two  azimuth  positions 
could  not  be  determined. 

The  variation  of  the  drag  moment  with  root  geometric  angle 
of  attack  is  presented  in  Figure  42.  It  can  be  seen  that  the  point  of 
minimum  drag  moment  is  roughly  the  same  for  both  azimuth  positions, 
approximately  6  to  7  degrees,  which  is  about  3  degrees  less  than  the 
angle  for  minimum  drag.  The  drag  moment  curves  are  again  approxi¬ 
mately  symmetrical  about  the  angle  for  minimum  moment  at  both  p  =  266 
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and  f  =  90  degrees.  The  rate  of  moment  increase  with  blAde  Angle  from 
the  minimum  is  significantly  larger  lor  p  Zbb  degree*  thAn  it  t*  lor 
p  -  90  degrees. 

Figure  43  presents  the  variation  of  the  ratio  ol  lilt  to  drag 
with  root  geometric  angle  of  attack.  It  i*  noted  that  the  maximum 
values  of  LID  occur  at  approximately  12  degree*  and  that  the  maxi* 
mum  value  at  p  -  266  degrees  is  much  1c**  than  It  i«  at  #  •  90  degree#. 

The  ratio  of  the  maximum  values  at  90  and  26t>  degree*  t*  approx!  - 
mately  3.  5:1,  which  compares  to  a  value  of  approximately  J.  b  I  baaed 
on  sectional  data.  The  fact  that  the  present  measurement  ol  thl*  ratio 
is  approximately  the  same  as  the  value  based  on  sectional  datA  l»  coin¬ 
cidental.  As  was  noted  previously,  the  individual  meA*ured  blade  force# 
did  not  exhibit  the  same  variations  with  angle  of  attack  a#  the  force#  ba#ed 
on  corresponding  two-dimensional  sectional  data. 

Figure  44  presents  the  variation  of  the  root  normal  force  And 
root  beamwise  moment  with  shaft  angle  for  an  azimuth  angle  ol  approxi¬ 
mately  180  degrees.  In  this  azimuthal  orientation,  the  blade  become# 
a  low-aspect:-ri-tio  wing  with  respect  to  the  free-slream  velocity.and 
the  shaft  angle  becomes  a  measure  of  the  geometric  angle  of  attack  lor 
the  wing.  As  a  result,  one  might  expect  that  the  variation  ol  the  blade 
aerodynamic  characteristics  with  shaft  angle  would  be  similar  to  that 
of  a  low-aspect-ratio  wing.  This  appears  to  be  the  case.A*  evidenced 
by  Figure  44.  It  is  noted  that  the  root  normal  force  I*  small  over  the 
entire  range  of  geometric  angle  of  attack.  Further,  the  variation  ol 
the  normal  force  with  geometric  angle  of  attack  is  nonlinrAr,  a*  might 
be  expected  for  a  low -aspect-ratio  wing.  The  corresponding  moment 
variation  is  also  nonlinear.  The  magnitudes  of  the  moments  reflect  the 
fact  that  the  moments  are  taken  about  the  blade  root  and  the  effective 
moment  arm  is  large.  The  moment  plot  indicates  that  blade  divergence 
might  be  reached  at  a  smaller  negative  shaft  angle  (i.  c.  ,  geometric 
angle  of  attack)  than  for  positive  shaft  angles.  It  is  believed  that  this 
is  due  to  the  effects  of  the  blade  spanwisc  flexibility.  At  zero  shaft 
angle  and  wind  off,  the  blade  tip  is  bent  down  approximately  |  foot 
due  to  the  blade  weight.  This  spanwisc  distribution  of  blade  deflection# 
causes  the  geometric  angle  of  attack  of  the  blade  with  respect  to  the 
free-stream  velocity  to  be  different  from  the  shaft  angle.  If  the  blAde 
were  rigid,  the  blade  geometric  angle  of  attack  and  thr  shaft  Angie  would 
be  identical.  Thus,  as  the  shaft  angle  is  increased  in  the  negative  direc¬ 
tion,  the  blade  flexibility  causes  a  more  rapid  increase  in  Angle  of  attack 
due  to  the  additive  contribution  of  the  blade  weight  and  aerodynamic  loading 
than  indicated  by  the  shaft  angle.  A  similar  effect  occur*  as  the  shaft 
angle  is  increased  in  the  positive  direction.  However,  at  some  positive 
shaft  angle,  the  aerodynamic  blade  load  and  moment  distribution  is  such 
as  to  just  balance  the  weight  load  and  moment  and  effectively  remove 
the  contribution  due  to  blade  deflection  due  to  weight.  Above  this  value 
of  shaft  angle,  the  blade  bending  due  to  aerodynamic  loading,  now  In  a 
positive  direction,  again  begins  to  contribute  in  an  additive  manner  to  the 
anf’le  of  attack. 


Thus,  if  divergence  were  to  occur,  it  might  be  expected  to 
occur  at  smaller  negative  shaft  angle  magnitudes  than  for  positive 
shaft  angles. 

Note,  also,  that  the  effects  of  blade  flexibility  tend  to  empha¬ 
size  the  nonlinear  behavior  of  the  blade  aerodynamics  with  geometric 
angle  of  attack  over  and  above  that  which  might  be  expected  if  the  blade 
were  rigid. 

Figure  45  presents  a  plot  of  the  spanwise  moment  distribution 
of  the  blade  for  various  shaft  angles.  Note  that  the  spanwise  distributions 
of  moments  for  negative  shaft  angles  are  very  similar,  while  those  for 
positive  shaft  angles  differ  in  shape  from  one  another.  The  shapes  of 
the  moment  distributions  for  positive  and  negative  shaft  angles  are  also 
different.  Not  until  acs  as  14  degrees  does  the  shape  obtained  for  positive 
shaft  angle  begin  to  look  like  the  shapes  obtained  for  negative  shaft  angles. 
It  is  believed  that  this  behavior  is  also  due  to  the  effects  of  blade  spanwise 
flexibility.  It  is  proposed  that  the  change  in  shape  of  the  moment  distri¬ 
bution  curves  in  the  range  -2 <aes*.  7  degrees  is  attributable  to  the  reduced 
effect  of  blade  spanwise  deformations  which  are  minimized  in  this  range 
due  to  the  balancing  of  the  aerodynamic  and  gravity  loads. 

The  fact  that  increased  dynamic  pressure  increases  the  aero¬ 
dynamic  loading  and,  hence,  the  blade  deflections  would  lead  to  the 
expectation  that  increasing  dynamic  pressure  would  increase  the  non¬ 
linear  behavior  of  the  blade  aerodynamics.  This  seems  to  be  substan¬ 
tiated  in  Figure  44  by  the  behavior  of  the  normal  force  and  moment  at 
?  =  45  psf  as  compared  with  their  behavior  at  q.  =  28  psf. 

Figure  46  presents  the  variations  of  the  root  edgewise  force 
and  moment  and  the  root  pitching  moment  with  shaft  angle.  It  is  noted 
that  all  of  these  quantities  are  not  much  larger  than  the  possible  measure¬ 
ment  errors  ;  hence,  there  is  considerable  scatter  in  the  data. 

Figures  47  through  51  present  the  azimuthal  variations  of  the 
root  normal  force,  root  beamwise  moment,  root  edgewise  force,  root 
edgewise  moment, and  root  pitching  moments,  respectively,  as  measured 
by  the  balance  system  for  three  different  shaft  angles  and  for  a  number 
of  different  blade  pitch  angles.  Steady  values  of  the  various  forces  and 
moments  could  not  be  obtained  for  certain  combinations  of  azimuth  and 
blade  pitch  angles  because  of  either  a  dynamic  instability  or  load  limits. 
These  areas  have  been  noted  on  all  the  figures.  The  limit  load  was 
based  upon  a  30,  000-psi  beamwise  stress  limit  that  was  arbitrarily 
imposed  at  the  critical  blade  section.  Because  of  the  beamwise  stress 
developed  by  the  gravity  load,  larger  positive  values  of  aerodynamic 
beamwise  moments  could  be  obtained  tnan  the  approximately  6,  000  ft -lb 
negative  aerodynamic  beamwise  moment  limit. 

No  data  are  presented  for  azimuth  angles  of  0  +  30  degrees 
because  of  interference  effects  introduced  by  the  wake  of  the  rotor 
head. 
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The  variation  of  the  root  normal  force  with  azimuth  angle  for 
zero  shaft  angle,  shown  in  Figure  47,  is  what  might  be  expected,  in 
that,  as  the  blade  angle  is  increased,  the  normal  force  on  the  advancing 
side  gets  larger  in  the  positive  direction  and  larger  values  in  the  nega¬ 
tive  direction  are  obtained  on  the  retreating  side.  It  can  also  be  seen 
that,  as  the  root  blade  angle  is  decreased,  negative  values  of  the  normal 
force  can  be  obtained  on  the  advancing  side  and  positive  ones  on  the 
retreating  side.  Regardless  of  the  blade  angle,  the  normal  force  is 
approximately  zero  at  the  180-degree  azimuth  position.  The  change 
in  the  azimuthal  variation  of  the  normal  force  that  is  obtained  with 
the  blade  pitch  angle  suggests  the  possibility  of  developing  a  constant 
normal  force  around  the  azimuth  if  the  proper  blade  pitch-control 
schedule  can  be  provided.  Since,  at  zero  shaft  angle,  the  normal 
force  is  zero  for  any  blade  pitch  angle  at  j V  -  180  degrees,  the  only 
possible  constant  normal  force  is  zero.  Figure  48  presents  the 
approximate  blade  pitch-control  schedule  that  would  be  required  at 
zero  shaft  angle  to  obtain  a  zero  normal  force  around  the  azimuth. 

This  curve  was  constructed  by  cross  plotting  the  data  presented  in 
Figure  47  to  determine  the  blade  angles  which  would  result  in  a  zero 
normal  force  at  various  azimuth  positions.  While  the  variation  pre¬ 
sented  in  Figure  48  must  be  considered  as  approximate,  these  results 
demonstrate  that  the  pitch-control  schedule  required  to  obtain  zero 
normal  force  at  the  blade  root  around  the  azimuth  would  be  rather 
complicated. 

The  results  presented  in  Figure  47  for  other  shaft  angles  show 
that  changes  in  the  shaft  angle  have  a  large  effect  on  the  azimuthal  varia¬ 
tion  of  the  normal  force.  The  exception  to  this  is  at  ^  =  180  degrees 
where  the  normal  force  has  only  a  small  variation  with  shaft  angle  and 
the  effect  of  blade  pitch  angle  is,  as  previously  noted,  nonexistent.  For 
a  positive  shaft  angle  of  5.  18  degrees,  it  appears  from  the  results  pre¬ 
sented  that,  in  the  range  of  azimuth  angle  30<^<180  degrees,  the  blade 
may  have  been  stalled  over  much  of  its  span  as  the  blade  angle  was 
increased  from  10.  19  to  13.  21  degrees.  The  normal  force  does  not 
increase  significantly  for  this  increase  in  blade  root  pitch  angle. 

The  variations  of  the  root  beamwise  moment  with  azimuth 
position  for  three  shaft  angles  and  several  blade  pitch  angles  are 
presented  in  Figure  49.  The  variation  of  the  root  beamwise  moments 
with  azimuth  angle  is  as  what  might  be  expected  on  the  basis  of  the 
variation  of  the  root  normal  force  for  the  same  test  conditions.  It  is 
noted  that,  at  y  =  180  degrees,  the  change  in  the  root  beamwise  moment 
with  shaft  angle  is  larger  for  negative  shaft  angles  than  it  is  for  positive 
shaft  angles,  although  the  change  in  the  normal  force  was  about  the  same. 
This  could  be  the  effect  of  blade  beamwise  bending, which  probably  shifts 
the  load  outboard.  It  is  apparent, when  the  variations  of  the  root  normal 
force  and  root  beamwise  moments  with  azimuth  angle  are  compared, that 
the  pitch  control  schedule  required  to  make  the  beamwise  moment  con¬ 
stant  around  the  azimuth  would  be  different  from  that  required  to  accom¬ 
plish  the  same  objective  with  the  normal  force. 
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The  variations  of  the  root  edgewise  force  with  azimuth  angle 
are  presented  in  Figure  50  for  three  different  shaft  angles  and  for  a 
number  of  different  blade  root  pitch  angles.  In  viewing  these  curves, 
it  must  be  remembered  that  this  force  is  oriented  parallel  to  the  blade 
root  chord.  Since  the  blade  root  chord  is  at  an  angle  of  attack  with 
respect  to  the  airstream,  the  root  edgewise  force  is  not  the  same  as 
the  drag  force  which  is  in  the  windstream  axis  system.  It  is  inter¬ 
esting  to  note  that,  generally,  the  root  edgewise  force  is  larger  on 
the  advancing  side,  30  <f  <180  degrees,  than  it  is  on  the  retreating 
side,  180«j*<330  degrees,  for  the  range  of  blade  angles  tested, and 
that  the  minimum  variation  in  the  force  with  azimuth  angle  was  obtained 
for  a  shaft  angle  of  -4.  77  degrees.  The  edgewise  force  is  not  generally 
in  the  direction  one  would  expect, in  that,  with  the  leading  edge  facing  the 
airstream,  the  edgewise  force  is  such  as  to  put  the  leading  edge  in  com¬ 
pression.  The  strain-gage  data  verify  this  result, as  the  edgewise  strain 
indicates  leading-edge  compression  along  the  entire  blade  span  for  posi¬ 
tive  pitch  angles  at  ^  =  120  degrees.  It  is  noted,  however,  that,  when 
the  root  normal  force  and  edgewise  force  are  resolved  to  obtain  lift  and 
drag,  positive  values  of  drag  (with  respect  to  the  airstream  reference 
system)  are  obtained  on  the  advancing  side  of  the  rotor  disc. 

The  variations  of  the  root  edgewise  moments  produced  by  the 
root  edgewise  forces  with  azimuth  angle  are  shown  in  Figure  51,  It  is 
noted  that,  as  with  the  root  edgewise  force,  the  variation  of  the  root 
edgewise  moment  is  generally  larger  on  the  advancing  side  of  the  disc 
than  it  is  on  the  retreating  side  and  that  the  minimum  variation  was  again 
obtained  for  all  blade  angles  at  a  shaft  angle  of  -4.  77  degrees. 

At  a  shaft  angle  of  0.  01  degree,  the  root  edgewise  moment  for 
for  a  root  pitch  angle  of  +  1.11  degrees  abruptly  reverses  in  direction 
between  =  100  and  130  degrees.  No  explanation  for  this  anomalous 
behavior  was  found. 

The  variation  of  the  root  pitching  moment  about  the  blade  quarter 
chord  with  azimuth  (Figure  52)  is  what  one  would  expect  --  small  over 
the  advancing  quadrant  of  the  disc  and  large  over  the  retreating  quad¬ 
rant  of  the  disc, where  the  moment  arm  is  relatively  large  because  of 
the  shift  in  the  center  of  pressure  from  the  quarter-chord  to  the  three- 
quarter-chord  position.  It  is  apparent  from  the  data  presented  that 
relatively  high  control  loads  can  be  expected  over  the  retreating  quad¬ 
rants  of  the  disc  as  a  rotor  system  is  being  slowed  to  a  stop  during  the 
conversion  cycle  from  a  helicopter  to  a  fixed-wing  aircraft. 

While  other  data  of  the  type  obtained  during  the  present  tests 
are  not  generally  available,  some  data  obtained  for  a  three -bladed  folding 
rotor  system  in  the  nonrotating  mode  are  available  for  comparison  (Refer¬ 
ence  9).  The  blades  that  were  tested  had  a  diameter  of  33  feet,  a  14-inch 
chord,  and  -9.  43  degrees  of  linear  blade  twist,  and  were  cantilevered  at  the 
root.  That  rotor  system  had  2-1/4  degrees  of  coning, whereas  the  present 
rotor  system  had  none.  Except  for  the  coning  of  the  rotor  system,  the 
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two  rotor  blades  had  similar  geometric  characteristics.  The  relative 
stiffness  of  the  blade  reported  in  Reference  however,  was  higher 
than  that  for  which  data  have  been  presented  herein. 

Figure  53  compares  the  azimuthal  variation  of  the  root  beam- 
wise  moment  as  measured  by  the  balance  system  during  the  present 
tests  with  that  obtained  by  means  of  strain  gages  at  the  root  of  the  blade 
that  was  reported  in  Reference  9.  It  is  noted  that  the  azimuthal  varia¬ 
tions  of  the  moments  are  very  similar.  It  is  somewhat  surprising,  how¬ 
ever,  that  larger  beamwise  bending  moments  were  obtained  with  the 
smaller  diameter,  smaller  chord  blade  of  Reference  9  than  with  the 
blade  used  to  obtain  the  data  reported  herein.  While  the  coning  angle 
would  be  expected  to  have  this  effect,  such  a  large  change  in  the  moment 
for  only  2-1/2  degrees  of  coning  would  not  be  expected. 

As  previously  noted  in  Figures  47  through  52,  data  could  rot 
be  obtained  for  various  combinations  of  blade  pitch  angles  and  shaft 
tilt  angles  in  the  azimuth  angle  range  225  *  l^/<300  degrees  because  of 
dynamic  instabilities. 

Figure  54  presents  basic  stability  boundaries  that  were  deter¬ 
mined  from  the  data  obtained  during  the  tests  at  various  shaft  angles. 

The  data  presented,  unless  otherwise  noted,  were  obtained  at  a  dynamic 
pressure  of  28.  95  psf, which  corresponds  to  a  Reynolds  number  of  approxi¬ 
mately  10°  per  foot.  As  indicated  in  Figure  54,  the  stable  region  was 
between  the  two  boundaries,  and  the  unstable  regions  were  above  the  upper 
boundaries  and  below  the  lower  ones. 

At  azimuth  angles  of  less  than  220  degrees,  instability  boundaries 
were  not  obtained  for  values  of  the  independent  parameters  within  rea¬ 
sonable  bounds.  At  ^  =  224  degrees,  blade  stress  limits  were  reached 
before  the  lower  instability  boundaries  c^uld  be  determined  for  shaft 
angles  <  2.  5  degrees.  When  the  instability  was  encountered,  a  low- 
frequency  blade  beam-bending  oscillation  could  be  seen  to  start  and 
increase  in  amplitude.  As  the  amplitude  of  the  low-frequency  oscilla¬ 
tion  built  up,  a  much  higher  frequency  blade  beam-bending  oscillation 
would  also  become  visible.  While  the  oscillatory  characteristics  were 
approximately  the  same  on  both  sets  of  boundaries,  the  low-frequency 
bending  amplitudes  would  become  larger  when  the  blade  was  at  negative 
angles  of  attack  (lower  set  of  boundaries)  than  when  the  blade  was  at 
positive  angles  of  attack  (upper  set  of  boundaries).  For  example,  the 
low-frequency  oscillatory  amplitudes  associated  with  the  lower  bound¬ 
aries  would  sometimes  be  allowed  to  reach  +  3  to  4  feet  by  visual  esti¬ 
mation  before  the  blade  angle  was  changed  to  stop  the  oscillatory  motion, 
while  +  2  feet  was  about  the  maximum  for  the  upper  boundary. 

The  solid  points  shown  in  Figure  54  denote  data  points  that 
were  obtained  for  £  =  43.  7  psf,  and  the  arrow  that  is  associated  with 
each  point  is  the  direction  in  which  the  blade  angle  was  being  changed 
when  the  instability  was  encountered.  The  symbol  shape  indicates 
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the  ^  =  28.  95  psf  curve  with  which  each  data  point  is  associated.  The 
point  at  -  224  degrees  and  <f.  -  43.  7  psf  was  obtained  at  a  shaft  angle 
which  was  different  from  any  of  the  cf.  -  28.  95  psf  data.  From  the  results 
that  are  plotted,  however,  it  can  be  seen  that  there  is  a  reasonably  large 
effect  of  p  ,  particularly  at  if>  -  224  degrees. 

Analysis  of  the  oscillograph  records  showed  that  the  higher- 
frequency  oscillations  involved  chordwise  and  torsional  motions  at  a 
frequency  approximately  equal  to  the  first  cantilever  torsion  mode  of 
the  blade  (Table  VIII),  and  the  low-frequency  oscillations  were  at  the 
first-coupled  flapwise  bending  mode  of  the  blade  (Table  IX). 

The  results  shown  in  Figure  54  indicate  that  shaft  tilt  angle, 
azimuth  angle,  and  blade  pitch  angle  were  all  independent  parameters. 
While  the  shaft  tilt  angle  is  a  significant  parameter,  it  is  noted  that  its 
effect  is  approximately  independent  of  the  azimuth  angle.  In  an  attempt 
to  determine  if,  in  fact,  the  shaft  angle  and  blade  pitch  angles  were 
independent  parameters  as  indicated  in  Figure  54,  the  total  angle  of 
the  blade  root  with  respect  to  the  airstream  direction  was  computed 
as  a  function  of  azimuth.  The  results  of  combining  these  angles  are 
shown  in  Figure  55.  As  can  be  seen,  all  the  data  points  for  each  sta¬ 
bility  boundary  at  a  given  azimuth  fall  within  a  band  of  approximately 
±  1  degree.  Since  all  the  boundaries  for  different  shaft  angles  become 
essentially  a  single  boundary,  the  shaft  angle  and  blade  root  pitch  angle 
are  not  independent  parameters,  and  the  primary  parameter  of  the 
instability  at  a  given  azimuth  location  is  just  the  root  geometric  angle 
of  attack  with  respect  to  the  airstream.  It  is  felt  that  the  scatter  at 
each  of  the  azimuth  locations  at  which  data  were  obtained  could  be 
reduced  somewhat  if  the  small  amounts  of  the  steady  elastic  torsional 
deflections  of  the  blade  were  taken  into  account.  This  has  not  been 
done,  as  it  is  believed  that  such  a  small  correction  would  not  change 
the  basic  conclusion  that  was  reached  from  the  results  presented  in 
Figure  55. 

While  the  combination  of  the  shaft  and  blade  pitch  angles  into 
an  equivalent  blade  root  angle  relative  to  the  airstream  indicated  that 
these  two  angles  were  not  independent  parameters,  the  plot  of  the  data 
indicated  that  the  composite  blade  angle  at  which  the  instability  occurred 
was  still  significantly  dependent  upon  the  azimuth  angle.  It  is  believed 
that  this  apparent  dependence  upon  azimuth  angle  is  associated  with  the 
bending  deflections  experienced  by  the  blade  in  the  vicinity  of  the  stability 
boundaries.  Unfortunately,  no  information  is  currently  available  which 
defines  the  effect  on  aeroelastic  stability  boundaries  of  lifting  surfaces 
when  large  loads  or  deformations  are  present.  However,  an  intuitive 
argument  (similar  to  the  one  presented  on  page  26)  for  the  case  in  hand 
would  imply  that,  if  large  blade  deflections  existed,  they  contributed 
to  the  effective  angle  of  attack  of  the  blade.  Further,  since  the  blade 
loadings  and,  hence,  the  blade  deflections  are  strong  functions  of  the 
azimuth  angle,  it  might  be  expected  that  the  instability  boundaries  would 
also  exhibit  a  similar  dependence. 
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Table  VI  IT 

UNCOUPLED  BLADE  CANTILEVER  FREQUENCIES 

MODE 

FREQUENCY 

cps 

MODE 

FREQUENCY 

cps 

1st  BEAM  BENDING 

1.18 

5th  BEAM  BENDING 

72.03 

2nd  BEAM  BENDING 

8.03 

1st  TORSION 

26.00 

3rd  BEAM  BENDING 

22.46 

2nd  TORSION 

69.27 

4th  BEAM  BEADING 

43.17 

3rd  TORSION 

113.56 

Table  JX 

COUPLED  BLADE  CANTILEVER  BENDING  FREQUENCIES 

MODE 

FREQUENCY 

cps 

GENERALIZED  COORDINATES  i  CHARACTERISTIC  SHAPE 

BEAM  BENDING 

SHAPE 

CHORD  BENDING 

SHAPE 

1st  COUPLED  BENDING 

1.23 

1.00 

1st  CANTILEVER 

-0.15 

1st  CANTILEVER 

2nd  COUPLED  BENDING 

7.08 

1.00 

2nd  CANTILEVER 

+  1.60 

1st  CANTILEVER 

3rd  COUPLED  BENDING 

8.76 

1.00 

2nd  CANTILEVER 

-0.62 

1st  CANTILEVER 

4th  COUPLED  BENDING 

23.22 

1.00 

3rd  CANTILEVER 

+0.58 

2nd  CANTILEVER 

5th  COUPLED  BENDING 

44.42 

1.00 

4th  CANTILEVER 

-0.18 

2nd  CANTILEVER 

6th  COUPLED  BENDING 

47.24 

1.00 

4th  CANTILEVER 

+5.10 

2nd  CANTILEVER 

7th  COUPLED  BENDING 

72.44 

1.00 

5th  CANTILEVER 

-0.04 

3rd  CANTILEVER 
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In  order  to  illustrate  the  sensitivity  of  the  boundary  to  blade 
pitch  angle,  Figure  58  presents  the  oscillatory  strain  data  for  the  blade 
at  a  pitch  angle  2  degrees  less  than  that  for  which  the  results  were  pre¬ 
sented  in  Figure  57.  As  can  be  seen  from  the  traces  presented  in  Fig¬ 
ure  57,  the  torsional  oscillations  are  fairly  persistent  while  the  chordwise¬ 
bending  traces  do  not  show  much  indication  of  oscillations  at  the  torsional 
frequency.  It  is  also  noted  that  there  is  not  much  low-frequency  beam 
bending  for  these  conditions.  A  prominent  characteristic  of  both  of  the 
instability  boundaries  was  that  high-frequency  oscillatory  motions  did 
not  become  apparent  visually  until  after  the  low-frequency  beamwise¬ 
bending  oscillations  became  fairly  large.  This  was  a  very  beneficial 
feature  during  the  tests,  as  the  stability  boundaries  were  deter  mined 
by  visual  observation  using  the  beamwise -bending  oscillations  as  an 
indicator. 

The  characteristic  frequency  of  the  instabilities  was  nearly 
equal  to  the  uncoupled  torsional  frequency  of  the  blade.  Also,  the 
local  geometric  angles  of  attack  of  the  blade  sections  were  large. 

These  factors  led  to  a  suspicion  that  the  instabilities  were  associated 
with  stall  flutter.  Evidence  that  supports  this  contention  is  presented 
below. 


Figure  56  presents  an  oscillograph  record  showing  the  beam- 
wise,  chordwise,  and  torsion  strain  signals  at  a  noted  instability  con¬ 
dition.  Although  the  traces  are  somewhat  weak  because  of  the  repro¬ 
duction,  the  following  characteristic  features  can  be  determined  from 
the  strain  traces: 

1.  The  torsional  oscillations  are  approximately  at 
the  uncoupled  torsion  frequency  of  the  blade 
(Table  VIII). 

2.  The  chordwise -bending  oscillations  are  at  the 
uncoupled  torsional  frequency. 

3.  Large  beamwise  oscillatory  bending  occurs  at 
the  first-coupled  bending  mode  frequency  of  the 
blade. 

4.  There  is  almost  no  high-frequency  response  at 
the  uncoupled  torsional  frequency  in  the  beam- 
wise  strain  signals. 

5.  The  high-frequency  torsional  and  chordwise 
oscillations  are  modulated  at  the  first -coupled 
flapwise -bending  frequency  (Table  IX). 
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Not  all  of  these  features  were  distinctive  in  all  of  the  oscil¬ 
lograph  records  taken  at  the  stability  boundary  points.  However,  the 
particular  records  discussed  herein  are  representative  of  those  in 
which  these  characteristics  prevailed. 

The  strain  signals  presented  in  Figure  56  are  for  a  point 
on  the  upper  stability  boundary.  For  this  point,  note  that,  in  addition 
to  the  above  features,  the  modulated  amplitudes  of  chordwise  bending 
and  torsion  reached  maximums  just  after  the  maximum  upswing  of 
the  beamwise -bending  signal  occurred,  and  the  chordwis e -bending 
oscillations  are  in  phase  with  the  torsional  oscillations. 

Figure  57  presents  the  oscillatory  strain  signals  for  a  point 
on  the  lower  instability  boundary  as  plotted  in  Figures  54  and  55.  It 
is  noted  that  the  same  basic  characteristics  as  listed  for  the  oscillatory 
strain  traces  shown  in  Figure  56  are  present,  except  that  the  chordwise¬ 
bending  signals  are  out  of  phase  with  the  torsion,  and  the  peak  ampli¬ 
tude  of  the  modulated  oscillation  occurs  just  after  the  maximum  down¬ 
swing  of  the  low-frequency  beam  benaing. 

To  relate  the  motions  of  the  blade  to  the  signals  presented  in 
Figures  56  and  57,  it  is  noted  that  the  positive  directions  of  beam  bending, 
chord  bending,  and  torsion  are  up  on  the  oscillograph  record.  Positive 
beam  bending  creates  compression  on  the  top  surface,  positive  chord 
bending  creates  compression  at  the  geome  ric  trailing  edge  of  the  blade, 
and  positive  torsion  twists  the  blade  in  a  nose-up  direction. 

Thus,  if  Figure  56  is  considered,  it  appears  that  the  blade 
beamwise -bending  oscillations  are  such  as  to  contribute  a  plunging 
velocity, which  increases  the  absolute  magnitude  of  the  local  angle  of 
attack  (remember  that  the  free-stream  direction  is  from  the  trailing 
edge  to  the  leading  edge)  while  the  blade  is  moving  up, and  decreases 
the  local  angle  of  attack  as  it  moves  down.  Thus,  if  the  blade  were 
close  to  the  stall  range  due  only  to  the  steady  angle  of  attack,  the 
beamwise -bending  oscillations  would  carry  the  blade  into  the  stall 
range  during  part  of  the  cycle  and  then  out  again  during  the  remainder 
of  the  cycle.  Thus,  the  bending  oscillations  could  precipitate  the  stall 
which  would  now  take  on  dynamic  stall  characteristics  rather  than 
static  stall  characteristics,  and  it  might  be  expected  that  the  torsional 
oscillations  would  grow  and  decay  in  phase  (nearly)  with  the  beamwise¬ 
bending  oscillation  frequency. 

On  the  lower  boundary  where  the  geometric  angles  of  attack 
are  in  the  opposite  direction,  it  would  be  expected  that  the  maximum 
amplitude  of  the  torsional  oscillations  would  occur  just  after  the 
maximum  downward  oscillation  of  the  beamwise  bending.  This 
appears  to  be  substantiated  in  Figure  57. 
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Other  experiments  precisely  like  those  reported  herein  have 
not  been  conducted  previously.  However,  there  have  been  a  number 
of  investigations  conducted  to  determine  regions  of  negative  pitch 
damping  associated  with  airfoils  oscillating  about  the  stall  angle  of 
attack.  These  can  be  referred  to,  at  least  in  a  qualitative  manner, 
to  determine  if,  for  the  conditions  of  the  present  tests,  one  might 
expect  a  stall-associated  instability.  It  should  be  noted,  however, 
that  for  the  present  tests,  the  airfoil  trailing  edge  was  into  the 
airstream,  and  the  aerodynamic  characteristics  of  the  airfoil  would 
be  expected  to  be  somewhat  different  from  those  for  an  airfoil  having 
a  normal  orientation  with  respect  to  the  airstream.  It  is  believed 
that  the  effect  of  the  "sharp  leading  edge"  would  be  a  difference  in 
the  static  aerodynamic  characteristics  around  stall,  and  the  difference 
in  the  dynamic  characteristics  of  stall  would  be  primarily  due  to  the 
effect  of  having  the  pitch  axi3  at  the  three-quarter  chord.  Halfman 
(Reference  10)  showed  that  the  primary  effect  of  leading -edge  sharp¬ 
ness  is  a  reduction  in  the  static  stall  angle.  Ham  (Reference  11) 
indicates  that  the  primary  effect  of  the  aft  pitch  axis  is  to  reduce  the 
dynamic  stall  angle.  For  the  reduced  frequency  of  0.  70  at  which  the 
instabilities  were  obtained  in  the  present  study,  the  data  that  Rainey 
presented  in  Reference  12  indicates  that  negative  aerodynamic  damping 
can  be  obtained  for  the  blade  angles  at  which  the  present  instabilities 
were  encountered.  This  fact  is  believed  to  further  substantiate  the 
thesis  that  the  present  instabilities  could  have  been  stall  flutter, 

E.  DEVELOPMENT  OF  A  LINEARIZED  AERODYNAMIC  PREDICTION 

TECHNIQUE - 


A  technique  of  predicting  the  aerodynamic  characteristics  of 
a  rotor  blade  at  arbitrary  azimuthal  position  has  been  developed  using 
conventional  lifting  surface  theory  tailored  for  this  specific  aoplication. 
The  detailed  development  of  the  theoretical  prediction  technique  and  the 
method  for  solving  the  resulting  equations  arc  presented  in  the  Appendix. 
The  wing  was  represented  by  a  distribution  of  bound  vorticity  and  the 
wake  by  trailing  vorticity,  whose  strength  distribution  is  related  to  the 
bound  vorticity  by  the  Helmholtz  laws.  Application  of  the  Biot-Savart 
law  to  calculate  the  velocity  field  associated  with  a  distribution  of 
vorticity  results  in  the  following  integral  equation  relating  the  local 
wing  loading  and  the  local  wing  slope  (i.  e.  ,  the  angle  between  the  local 
tangent  to  the  mean  camber  line  and  the  wind  reference  in  the  plane  of 
the  wind  vector). 
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where  ofx.y)  is  the  local  wing  slope  and  is  the  local  wing  loading 

(difference  between  upper  and  lower  surface  pressures).  Here  the  Kutta- 
Joukowski  relation  relating  loading  to  the  bound  vorticity  has  been  used, 
i.  c  *it)'pVY(x..y.)  where  y(x..<j.)  *»  «he  bound  vorticity  distri¬ 

bution.  The  derivation  of  the  integral  equation  (Equation  I)  is  given 
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in  References  3  and  4.  In  this  equation,  ac(x,y)  is  the  known  quantity 
and  is  the  unknown.  The  method  of  solution  employed  is 

that  of  assumed  loading  modes.  That  is,  the  functional  form  of  £(x0,yc) 
is  assumed,  this  form  being  linear  in  certain  undetermined  constants. 
These  constants  are  later  evaluated  by  requiring  exact  tangency  at  a 
specified  number  of  points  on  the  wing.  However,  before  this  was 
accomplished,  it  was  advantageous  to  perform  a  transformation  of 
variables  on  Equation  (1).  This  was  done  by  choosing  a  coordinate 
system  fixed  to  the  blade  with  origin  at  the  geometric  center  of  the 
blade.  The  spanwise  dimension  was  nondimensionaiized  by  onc-half 
of  the  span  and  the  chordwise  dimension  by  onc-half  of  the  chord.  This 
transformation  deforms  the  wing  to  a  yawed  square  and  is  depicted  in 
the  following  sketch. 


Equation  (I)  is  then  transformed  to 
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where  *t0'(g0.  ?*)i3  assumed  to  be  of  the  following  form: 
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where  thr  di  and  V^/are  the  undetermined  constants  to  be  evaluated.  It 
is  noted  that  the  first  bracketed  expression  in  and  fk  gives  the  type 
of  loading  that  would  be  expected  if  one  could  formally  split  the  span* 
wise  and  chordwise  flow  over  the  wing.  Thece  were  chosen  in  this 
form  so  that  the  loading  would  approach  the  proper  limits  for  f  equal 
to  aero  or  90  degrees.  The  last  bracketed  term  in  f,  and  ta  was  added 
to  prevent  the  loading  at  the  apex  from  dropping  to  aero  for  intermediate 
values  of  f».  It  will  be  noted  that  the  loadings  were  formed  by  multi¬ 
plying  the  first  thr«*e  trrms  of  a  Dirnbaum  chordwise  series  by  an 
Appropriate  spanwise  function.  The  composite  loading  then  possesses 
an  appropriate  airfoil-type  singularity  at  the  leading  edges  And  has  no 
loading  at  the  trailing  edges  to  conform  to  the  Kutta  condition.  The 
expression  for  the  loading  contains  36  undetermined  constants. 
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The  expression  for  J?a  is  then  substituted  in  Equation  (2),  which 
can  then  be  integrated  to  give  an  expression  of  the  form 

06  (*•»!)  ’  *  *-(*■'?') (6) 

which  is  linear  in  the  a.’s  and The  cl’s  and  are  now  found  by 
selecting  j  values  of  /  and  i  values  of  y  such  that  ixj  -  36.  In  this 
fashion,  a  system  of  36  linear  algebraic  equations  is  obtained  from 
Equation  (6).  This  system  of  equations  is  then  solved  for  the  as  and 
J,'s .  The  integration  required  to  obtain  Equation  (6)  must  be  performed 
numerically.  Since  the  kernel  of  Equation  (2)  contains  a  strong  singu¬ 
larity,  this  involves  considerable  work  in  isolating  the  singularities 
before  numerical  integration  can  be  performed.  The  integration  is 
considerably  eased  by  making  use  of  Chebyshev-Gauss  quadrature. 

The  theory  for  this  quadrature  scheme  is  given  in  Reference  5.  The 
application  of  this  quadrature  to  lifting  surface  problems  was  evidently 
first  made  by  Hsu  (Reference  6).  To  a  large  extent,  the  present  analysis 
is  an  adaptation  of  Hsu's  technique  to  the  yawed  wing  case.  However, 
since  the  wing  planform  for  the  present  application  was  not  symmetrical 
for  arbitrary  yaw,  many  of  the  computational  economics  available  for 
the  symmetrical  case  could  not  be  employed. 

The  loading  modes  used  in  the  present  analysis,  as  can  be  seen 
from  Equations  (4)  and  (5),  were  essentially  parabolic  in  the  spanwise 
direction.  The  conventional  choice  for  the  functional  variation  in  the 
spanwise  direction  is  elliptic.  However,  this  choice  is  excluded  from 
the  present  analysis  because  certain  integrations  to  be  performed  require 
a  finite  spanwise  slope  of  the  wing  tip. 

For  symmetrical  planforms  and  elliptic  spanwise  mode  shape, 
Hsu  (Reference  6)  was  able  to  use  the  theory  of  Gaussian  quadrature 
to  develop  explicit  expressions  for  selection  of  the  spanwise  and  chord- 
wise  collocation  points  such  that  the  error  in  lift  is  minimized  for  a 
given  number  of  points.  Unfortunately,  the  same  analysis  cannot  be 
employed  for  the  yawed  case. 

Since  the  combined  loading  on  the  wing  is  of  the  form  given  by 
Equation  (3),  it  was  thought  to  be  desirable  to  use  the  same  expression 
to  obtain  the  chordwise  and  spanwise  collocation  points  and  that  these 
points  should  be  symmetrically  located  about  the  midchord  and  midspan 
line  of  the  wing.  The  selections  were 
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Location  of  the  collocation  points  at  these  chordwise  coordinates  results 
in  a  minimum  error  (for  a  given  number  of  chordwise  collocation  points) 
in  the  chordwise  integration  for  sectional  lift.  The  selected  spanwise 
collocation-point  locations  provide  the  minimum  error  in  the  deter¬ 
mination  of  the  moment  of  a  parabolic  spanwise  loading. 

Since  the  present  analysis  involves  several  departures  from 
conventional  lifting  surface  developments,  it  was  felt  desirable  to 
analyze  the  zero  sweep  or  yaw  case  using  the  same  techniques  used 
in  the  yawed  case  and  to  compare  the  results  with  References  6  and  7. 

In  particular,  it  was  desired  to  assess  the  influence  of  using  parabolic 
rather  than  elliptic  spanwise  loading  modes  and  alternative  selections 
of  collocation  points.  This  has  been  done  for  an  aspect-ratio-one, 
rectangular  wing  and  the  results  are  summarized  in  Table  X.  The 
present  program  duplicates  Hsu's  results  (Reference  6)  when  elliptic 
spanwise  loading  modes  are  used  with  Hsu's  choice  for  the  locations  of 
collocation  points.  With  the  symmetrical  distribution  of  Equation  (7) 
as  the  choice  of  collocation-point  locations,  the  elliptic  loading  modes 
give  results  which  compare  very  favorably  with  those  of  References  6 
and  7.  With  the  same  choice  of  locations  for  the  collocation  points,  the 
parabolic  spanwise  loading  modes  give  a  total  lift  that  is  only  slightly 
high.  However,  the  resulting  spanwise  distribution  of  integrated  chord- 
wise  loading,  in  this  case,  peaks  near  the  wing  tips.  Therefore,  the 
representation  using  parabolic  loading  modes  in  combination  with  15 
collocation  points  located  according  to  Equation  (7)  is  considered  unsat¬ 
isfactory.  In  an  attempt  to  improve  the  results  with  parabolic  loading 
modes,  a  minimum-square-error  procedure  was  used  with  30  collo¬ 
cation  points  and  15  undetermined  constants.  This  did  not  result  in 
significant  improvement.  It  appears  then  that  more  modes  and/or  an 
alternate  location  scheme  for  the  collocation  points  will  have  to  be  con¬ 
sidered  to  achieve  a  satisfactory  representation  employing  parabolic 
spanwise  loading  modes. 

Although  there  was  some  question  of  the  adequacy  of  the  para¬ 
bolic  loading  modes  and  the  collocation-point  locations  in  the  unyawed 
case,  the  computer  program  for  the  arbitrarily  yawed  wing  was,  never¬ 
theless,  developed  using  this  representation.  However,  its  application 
to  the  case  of  an  aspect-ratio-one  wing  yawed  at  45  degrees  did  not 
predict  the  symmetries  in  the  loading  which  are  inherent  in  the  problem. 
This  failure  cannot  be  attributed  to  the  use  of  parabolic  modes  or  par¬ 
ticular  collocation -point  locations,  but  whether  it  is  due  to  an  error  in 
the  theoretical  formulation  or  in  the  computer  program  remains  an 
unanswered  question  until  further  study  can  be  made. 

F.  CORRELATION  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


In  view  of  the  difficulties  discussed  above,  theoretical  calcu¬ 
lations  for  comparison  with  measurements  were  limited  to  the  special 
cases  for  the  unyawed  wing  corresponding  to  the  rotor  blade  at  f  =  90, 
180,  and  270  degrees.  Elliptic  spanwise  loading  modes  were  used  with 


39 


Table  I 

COMPARISON  OF  VARIOUS  THEORETICAL  RESULTS  FOR 
LIFT  ON  AN  ASPECT  RATIO  =  1.0  RECTANGULAR  W I NG 


CAL  ELLIPTIC  LOADING  1.4535 
MODES  M,,  =  0 


(2J-l)fT 

1 1  =  COS  —T. - 5  Jmax  =  3 


HSU  (REF.  6)  Mn  =  0.2  1.51 

^  =  0  1.48 


(21—1 )  rr 

\  =  COS  — T~.  ’  'max  “  ® 

‘‘max 


r  2Jrr  .  i  -  o 

4 1  =-COS -  ’  Jmax  ~  6 

3-imax+l 


HSU  (REF.  6)  Mn  =  0.2  1.497 

Mu  =  0  1 . 47 


*i  =-c°S  - - Tp  'max  =  3 

'max  1 


r  2jrr 

S,  =-C0S  - - —  ;  jmav  =  5 

J  2J  +i  max 

4Jmax  1 


7  j  =-C0S 


’  'max  =  ® 


NASA  KERNEL  FUNCTION  1.455 
(REF.  7) 


(FOR  SELECTION  OF  COLLOCATION 
POINTS  SEE  REF.  7) 


1.57  SLENDER-WING  THEORY 
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18  collocation  points  located  in  accordance  with  Equation  (7).  The 
computed  results  and  the  corresponding  measured  values  are  pre¬ 
sented  in  Table  XI  for  the  three  azimuth  positions  when  the  shaft  tilt 
angle,  <Xa  ,  =  0  degrees,  the  blade  root  pitch  angle,  ,  =  +  13  degrees, 
and  the  dynamic  pressure,  ^  ,  =  28,  85  psf.  In  general,  the  agreement 
between  the  theoretical  and  experimental  results  is  not  satisfactory. 
However,  as  discussed  below,  some  of  the  discrepancies  can  be 
explained. 


Table  XL 

COMPARISON  OF  MEASURED  AND  PREDICTED  RESULTS 

ccs  =  0°,  0^  -  +43°,  <}  =  26.85  psf 

■f  -  90  DEGREES 

f  =  180 

DEGREES 

y-=  270  DEGREES 

MEASURED 

COMPUTED 

MEASURED 

COMPUTED 

MEASURED 

COMPUTED 

LIFT 

(lb) 

700 

520.11 

3 

5.44 

-333 

-529.2 

ROLLING  MOMENT 

6074 

4246.7 

-61 

-104.74 

-3536 

-4402.5 

(ft- 1b) 

PITCHING  MOMENT 

91 

-2.705 

-1 

-1.651 

238 

339.26 

(ft- lb) 

DRAG 

22 

11.44 

0 

-34 

-13.17 

(lb) 

DRAG  MOMENT 

331 

64.88 

0 

-518 

-25.57 

(ft- lb) 

In  comparing  the  results  obtained  at  ^  =  90  degrees,  it  is 
noted  that  the  lift,  rolling  moment,  drag  force,  and  drag  moment  are 
underpredicted.  It  is  believed  that  a  reason  that  the  drag  force  and 
moment  are  less  than  measured  is  the  fact  that  the  profile  drag  and 
friction  drag  were  not  considered  in  the  calculations. 

The  difference  in  the  predicted  and  measured  lift  and  rolling 
moment  corresponds  to  the  incremental  lift  and  rolling  moment  that 
would  be  obtained  from  a  2-1 /2-degree  change  in  the  angle  of  attack. 
As  previously  noted  in  this  report,  the  experimental  results  that  were 
obtained  indicated  that  the  root  geometric  angle  of  attack  for  zero  lift 
was  approximately  1-1/2  degrees  different  from  that  expected  for  a 
rigid  blade  having  no  tip  effects.  Flow  angularity,  however,  is  not 
believed  to  be  the  cause  of  the  discrepancy  between  theory  and  experi¬ 
ment,  as  a  flow  angularity  of  anywhere  near  the  magnitude  required 
could  not  be  justified  on  the  basis  of  calibration  results  obtained  for 
the  NASA -Ames  40-by-80-foot  wind  tunnel. 

For  f  =  270  degrees,  the  same  general  comments  concerning 
the  correlation  of  the  theoretical  and  experimental  results  apply. 
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CONCLUSIONS 


It  is  concluded  that  the  primary  objective  of  the  research  pro¬ 
gram  was  realized,  i.i  that  a  set  of  blade  strain  and  aerodynamic  loading 
data  were  obtained  simultaneously  for  a  nonrotating  helicopter  blade 
arbitrarily  oriented  with  respect  to  the  free  stream.  In  addition  to 
this  general  conclusion,  the  following  specific  observations  are  made 
as  regards  the  aerodynamic  characteristics  of  a  rotor  blade  fixed  with 
respect  to  the  free  stream: 

1.  At  azimuth  angles  of  approximately  90  and  270 
degrees,  the  variations  of  the  lift,  rolling  moment, 
and  pitching  moment  about  the  quarter  chord  with 
angle  of  attack  are  as  expected,  based  on  sectional 
data. 

2.  The  minimum  drag  forces  at  azimuth  angles  of  90 
and  270  degrees  are  about  equal,  which  is  not  what 
would  be  expected,  based  on  sectional  data. 

3.  The  normal  fo^co  and  rolling  moment  tend  to  zero 
at  an  azimuth  angle  of  180  degrees  for  all  shaft 
tilt  and  blade  angles  tested. 

4.  The  data  appear  to  be  consistent  with  the  expectation 
that  blade  beamwise  flexibility  would  reduce  either 
the  dynamic  pressure  or  the  rotor  shaft  tilt  angle 

at  which  blade  bending  divergence  occurs  at  an 
azimuth  angle  of  180  degrees. 

5.  The  blade  pitch  angle  is  an  ineffective  control 
parameter  at  an  azimuth  angle  of  180  degrees. 

6.  For  all  azimuth  angles  except  those  at  180  +  10 
degrees  and  360  t  10  degrees,  the  rotor  shaft 
angle  and  blade  pitch  angle  are  significant  control 
parameters. 

7.  The  azimuthal  variations  of  the  drag  and  drag  moment 
for  any  blade  pitch  setting  decrease  as  the  shaft  tilt 
angle  is  varied  from  +  5  degrees  to  -  5  degrees. 

8.  The  root  normal  force  or  root  beamwise  bending 
moment  can  be  made  zero  at  each  azimuthal 
position  by  proper  choice  of  blade  angle. 

9.  The  magnitude  and  azimuthal  variation  of  the  blade 
root  pitching  moment  about  the  blade  root  quarter 
chord  tend  to  be  small  for  the  blade  on  the  advancing 
side  but  large  for  the  retreating  side. 
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10.  In  the  azimuth  angle  range  of  225  S  <  330  degrees, 
a  dynamic  instability  is  encountered  for  both 
positive  and  negative  blade  root  geometric  angles. 
There  is  some  evidence  in  the  records  presented 
that  the  instabilities  have  the  characteristics  of 
stall  flutter. 

11.  The  contribution  of  the  flexible  blade  bending  slope 
to  the  effective  angle  of  attack  is  significant  in  the 
establishment  of  the  instability  boundary. 


RECOMMENDATIONS 


On  the  basis  of  the  results  obtained  during  this  program,  it 
is  recommended  that  the  following  additional  research  efforts  be  under¬ 
taken  in  support  of  the  successful  development  of  stoppable  rotor  con¬ 
figurations  : 

1.  Obtain  aerodynamic  and  strain  data  for  a  stoppable 
rotor  configuration  in  the  stopping  and  starting 
operational  modes. 

2.  Conduct  additional  tests  with  rotor  blades  having 
different  stiffnesses,  twist  distributions,  and 
planforms  to  determine  the  effects  of  these 
parameters  on  the  instabilities  that  were 
encountered. 

3.  Continue  the  development  of  a  suitable  theory 
for  predicting  the  aerodynamic  characteristics 
of  stoppable  rotor  configurations. 

4.  Develop  a  theory  for  predicting  the  stall  flutter 
characteristics  of  stoppable  rotor  configurations. 
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Figure  1.  TEST  INSTALLATION  IN  WIND  TUNNEL 


Figure  2.  DIMENSIONAL  DRAWING  OF  TEST  INSTALLATION 
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Figure  3.  INSTALLATION  SHOWING  ACTUATOR  DRIVE  SYSTEMS 


Figure  4.  BALANCE  ASSEMBLY  ON  SUPPORT  STRUCTURE 
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Figure  5.  DIMENSIONAL  DRAWING  OF  FORCE  BALANCE  SYSTBI 


Figure  6.  GEOMETRIC  PROPERTIES  OF  THE  TEST  BLADE 


WEIGHT  PER  INCH,  lb/in. 


Figure 


Figure  9.  BLADE  BENDING  STIFFNESS  DISTRIBUTION 
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Figure  12.  PRESSURE  TAP  AND  STRAIN  GAGE  INSTALLATION  AT  85-PERCENT  SPAN 
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Figure  14.  AZIMUTH  AND  SHAFT  ANGLE  TRANSDUCERS 
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Figure  15.  INSTRUMENTATION  BLOCK  DIAGRAM 
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Figure  16.  INSTRUMENTATION  SETUP  IN  CONTROL  ROOM 
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Figure  17.  DATEX  II  STRIP  CHART  READOUT 
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Figure  18.  OPERATOR'S  CONTROL  PANEL 
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Figure  21.  SHAFT  TILT  ANGLE  CALIBRATION  CURVE 


INSTRUMENTED  BLADE  PITCH  ANGLE,  (9.,  deg 


Figure  22.  PITCH  ANGLE  CALIBRATION  CURVE  -  INSTRUMENTED  BLADE 
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Figure  24.  PRESSURE  GAGE  CALIBRATION  EQUIPMENT 


PRESSURE,  psi 


AMPLIFIER  SIGNAL,  volts  x  10 


Figure  25.  TYPICAL  CALIBRATION  CURVE  FOR  1-PSI  DIFFERENTIAL  GAGE 
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Figure  27.  TYPICAL  CALIBRATION  CURVE  FOR  4-PSI  DIFFERENTIAL  GAGE 
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28.  TYPICAL  CALIBRATION  CURVE  FOR  8-PSI 
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GAGE  CALIBRATION  CONSTANTS 


2.1 


TORSIONAL  STRAIN  GAGE  CALIBRATION  CONSTANTS 


BEAMWI SE  DEFLECTION,  in 
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Figure  34.  BLADE  BEAMWISE  DEFLECTION  UNDER  ITS  OWN  WEIGHT 
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EDGEWISE  DEFLECTION,  in 
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Figure  35.  BLADE  EDGEWISE  DEFLECTION  UNDER  ITS  OWN  WEIGHT 
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Figure  36.  TEST  OPERATOR'S  STATION  IN  WIND  TUNNEL  CONTROL  ROOM 
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Figure  37.  TYPICAL  PRINTOUT  FOR  EACH  DATA  POINT 
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Figure  38.  BLADE  LIFT  VS  ROOT  GEOMETRIC  ANGLE  OF  ATTACK  AT  p  =  90  AND  266  DEGREES 


BLADE  ROLLING  MOMENT  VS  ROOT  GEOMETRIC  ANGLE  OF  ATTACK 
VS  f  =  90  AND  266  DEGREES 
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BLADE  PITCHING  MOMENT  VS  ROOT  GEOMETRIC  ANGLE  OF  ATTACK 
FOR  #  =  90  AND  266  DEGREES 
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Figure  43.  BLADE  LIFT  FORCE/ BLADE  DRAG  FORCE  VS  ROOT  GEOMETRIC 
ANGLE  OF  ATTACK  FOR  f  =  90  AND  266  DEGREES 
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Figure  45.  BENDING  MOMENT  DISTRIBUTIONS  AT  ^  =  180  DEGREES  FOR  VARIOU 


ROOT  PITCHING  MOMENT,  ft-lb  ROOT  EDGEWISE  MOMENT,  ft-lb  ROOT  EDGEWISE  FORCE,  lb 


Figure  48.  VARIATION  OF  BLADE  ROOT  PITCH  ANGLE  FOR  ZERO  NORMAL  FORCE 
AT  ZERO  SHAFT  ANGLE 


94 


EDGE VI  SE  MONEVT 


■  ft. If 


INSTABILITY  AREA 


AZIMUTH  ANGLE,  /,  deg 


Figure  51.  ROOT  EDGEWISE  MOMENT  VS  AZIMUTH  ANGLE  FOR 
VARIOUS  SHAFT  ANGLES  AT  q  =  28.93 
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APPENDIX 

DETAILED  DEVELOPMENT  OF  THE  LIFTING  .SURFACE 

INTEGRAL  EQUATION  AND  ITS  SOLUTION - 


The  basic  lifting-surface  problem  is  embodied  in  the  solution 
to  the  following  integral  equations  for  £(%0,yo),  the  blade  loading: 


cz 


(*-*o) 


[(x-Z^f  +  ({/-%) 


‘] 


dxe  d,tjc 


(8) 


This  equation  is  derived  in  References  3  and  4.  The  "finite  part"  of  the 
divergent  integrals  is  understood  here  (see  References  3  and  4).  a.(z.y) 
is  the  local  slope  of  the  wing  surface  in  the  streamwise  direction  and  is 
considered  as  known.  S  is  the  wing  planform.  Equation  (8)  is  the  result 
of  representing  the  rotor  blade  by  a  distribution  of  bound  vorticity  [whose 
strength  is  proportional  to  l?(z0,e/e)  ]  and  the  wake  by  trailing  vorticity, 
whose  strength  is  related  to  the  bound  vorticity  by  the  Helmholtz  laws. 
Equation  (8)  can  then  be  interpreted  as  a  linearized  version  of  flow 
tangency  at  the  blade  surface.  Equation  (8)  is  solved  by  assuming  a 
functional  form  of  (%<,-&), which  is  linear  in  certain  unknown  constants. 
The  integration  indicated  is  then  performed  to  obtain  a  linear  relation 
between  y )  and  these  unknown  constants.  The  unknown  constants 
are  then  evaluated  by  requiring  exact  flow  tangency  at  a  compatible 
number  of  points  (referred  to  as  collocation  points)  in  the  rotor  plan- 
form.  The  integrations  must  be  performed  numerically  and  require 
isolation  of  the  strong  singularity  in  the  kernel  of  Equation  (8). 


The  overall  approach  to  the  problem  was  that  employed  by 
Hsu  (Reference  5).  However,  the  details  of  the  problem  are  vastly 
complicated  owing  to  the  nonsymmetrical  planform  under  consideration. 
Because  of  the  rectangular  planforms  that  are  under  consideration  at 
arbitrary  yaw  angles,  it  is  advantageous  to  make  the  following  trans¬ 
formations  of  independent  variables,  the  first  being  a  rotation  of  axis 
such  that 

/* 

X  -  £  cos  ft  -  y  'si* 1  ft 
* 

y  =  $  ft  -  y  cos  ft 


then,  a  nonuniform  stretching  by  letting 


I 


I 


where  c  i3  one-half  of  the  chord  of  the  wing  and  6  is  one-half  of  the  wing 
span.  As  can  be  seen,  this  transformation  reduces  the  wing  to  a  yawed 
square.  Equation  (8)  then  becomes 


,  be  /' /'  _  f  ctos  1. 

~3rrJ_f  [csMitfr-jr)+/>cosit'(7-f,)]£  /  IcVf-fJ+bty-fof]*  i 

It  is  now  assumed  that  is  of  the  following  form  : 


P,  Pi 

‘JR?  7^? 


(10) 


where 

Pi  *  ( fk,  +  0*0/  rjo  *  &oj  &0+  *? o  \  (  t  ~  4e) 

+  [al0  +  <z1t  rf0  +a,2%  +  a, j  ?*+  a,(4  7/  ]  ( /  -  {* ) 

[k<?  4  aH^o  4  O'  /Z  +  &ZS  *?o  4  &Z4  'll 


Z  *f  1°  J  G'SO  )  4  (f~€o  )  +  0'3z(t‘(o  )f»\  (11] 


Pz  *(l~£o  )  4-^0i  %0  4^t 9+€ff\(  t  ~  fy) 


4  [J»  +AzZot+Ju  *S+J,+fS  ]</’*/) 

4  [-^za  +  'tzt  £0  *  •J’/t  3o  4-^t3  £4  4^r+  (o  J  ( f'  1o  )  >?c  | 

*(i  “!*•  ~+f*)^*o(/-'?o)4Jjt('-'?S)  *  (12) 

This  form  then  contains  18  *'■*  and  18-^i,  for  a  total  of  36  unknown 
coefficients  for  which  a  solution  must  be  obtained.  These  mode  shapes 
also  restrict  the  range  of  V'  such  that  7  =  -1  and  (  =  -1  must  always 
physically  represent  the  leading  edges  of  th"  wing, since  these  lines 
are  the  square  root  singularity  lines  of  the  loading  modes.  This 
implies  that  the  azimuth  angle  must  be  in  the  range  of  Oipz-Z, 

However,  this  represents  no  loss  of  generality,  since  any  yaw  position 
can  be  transformed  to  an  equivalent  aerodynamic  problem  with  p  in 
this  range. 


I 
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Equations  (10),  (11), and  (12)  can  now  be  substituted  in  Equa 
tion  (9)  and  integrated  to  yield  an  expression  of  the  form 

&(£,  y  )  -  AOQ  &oo  ^  A  01  a,#/  +  A  02  &02  +  A  00  &oj  -f-  A  0+  ds0+  +  A  to  O-fo 
■+  A,,  a-,,  +  A u  a-u  A,3  cii3  +  A/4  a. 14  4  Azoato  +  A 2,0-2, 


4  A 22  a'2z 

+  A 20  0-23 

■*"  A 00  J00 

4  A0,^a, 

+■  A  if  Jit, 

4  L,2  Jf2 

+  A 22  A 22 

+  A  23  J2J 

“*■  A 24  CL24  A  jo  a. jo 

-*■  l-0zJ>02  -t-  A-0J  J-o3 

^  A, 3  Ji3  *■  Lf+1,4  + 

4  ijolje  + 


Aj/O-j/^  Aj 2  A32 

L04-loq-+  L-fO't'IO 

t-ZO'ttO  +  A  2/  -J-2) 
Aj/Jj,  +  L32/jz 


+  Ajj  Jjj  (13) 

where 

Axx  “  AAX(<r,  y) 

A- xx  *  Ajx  ($, 

Since  <x(f/j)  is  known,  j  values  of  {  and  i  values  of  7  such  that  t *j  =  36 
can  be  selected  and  exact  tangency  required  at  these  locations  (collocation 
points).  Then,  Equation  (13)  represents  a  system  of  36  linear  algebraic 
equations  that  can  be  solved  for  the  a/s  and  Us  . 

To  obtain  the  As  and  As  ,  the  integration  of  Equation  (9)  is 
divided  into  two  parts.  We  first  let 


q  !  ^  C  COS  -£,)  -  b  sin 

and 

?o  “  **  *  c  Co*t  ('?''?•) 

fo  "  f  *-§  ) 

Then,  c&n  be  expressed  as 

&((,  r?)  *  I,  +It 


(M) 

(15) 

(16) 
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where 


and 


T'  *  ■«?  Ut 

/  />  rV  P*fc-?s)  G-/tfe.n*>f-n)  , 

:*M~8rr  C  ll  {7^7  [io  -  So]  Z  s<o  > 


(17) 


(18) 


Only  the  detailed  evaluation  of  X  will  be  discussed  in  detail  since  there 
are  many  similarities  involved  in  the  two  integrations.  Pt  and  G,  are 
regular  functions  in  and  present  no  problems  for  integration. 

However,  the  square  singularity  in  the  denominator  of  the  kernel  func¬ 
tion  requires  special  treatment  before  numerical  calculation  can  be 
accomplished.  X,  is  first  integrated  with  respect  to  rj0 . 


/  c  r'f*  PtU.-'?*)  G,(<o‘ %•{.'?)  . 

'  ‘  -gr  Sll  JT-gW-bY  7‘  ** 

t  c  7*)  ~Pi((o,r?o)  G,(to  •  Ho  *  H)  . 

-fR?  [7j0  -  7oy  cos*+  I'*!0 

'  c  r’  4 «*•?•)  r' g, (f. .  .. 

'drr  b  /,  77-  //  J,  C«V  (19) 


b&  r  ^  (f  -?oyoost[{l-lts)+$  (r?- O'Y* 

<, -  b  Cos  +  [q,  -rf,]  ^  {C-(0)cosi//\j<.n £  cosfiq*  t)] 

( f-fo )  *  cos? 

1 t’fo)  COS  f/  \sui  'H(  So)*-ZCOi'l'(fl* 

The  (0  integration  of  the  last  integral  in  Equation  (19)  is  accompliahcd 
numerically  by  Chebyshev-Gauss  quadrature  (Reference  5).  Briefly, 
this  quadrature  can  be  summarized  as 


3 


Sf  ( f*  >  f.  q ) 


(20) 


where  the  xf  arc  selected  according  to  the  rule 


(21) 


»  • 
*4 


(22) 


108 


whore  the  integral  is  exact  if  f(x)  io  a  polynomial  of  degree  less  than 
2n-/.  Otherwise,  the  degree  of  precision  is  2*-f. 


Then 


8r  h  J. 


1  (fo  >  Ho  ) 


f  tf)  dfe 


•  -j-Z  rUp.i.)  ]/{*.(.  i) 

rr  P.t 


(23) 


and  it  is  to  be  remembered  that 


7*  *  ?•  ((*•  (•  n) 


Now,  continuing  treatment  of  the  first  integral  in  Equation  (19),  we  let 


P,((*  .  7#  )  -  n. ) 

V,) 


?•) 


(24) 


and, using  the  first  term  of  its  Taylor  Scries  expansion,  P,(?0.'}0  it 
calculated  as 


o) 


PP, 


n»a7» 


(25) 


Then Hi(f0.^0)  is  a  well-defined  function.  Equation  (19)  can  now  be 
written  as 


-  L  c  [' f*  ft)  &»((*•  ?*•  (•  7)  j  ,,  /  r  f*t.  -  ,  , 

Ir  frr  >4,1,  f)) 


t  c  r  r  1* )  i  o •  7*  )J  f'f®.  7*'  ^  7) 

•1  .( 


8rr  b 


cos'+ii7*;  (%•%) 


d%  d(. 


1  c  /'  *,(€.. 9.)  /'  C,(f0. _ 

•r 4  J,  /'-<•/  ,/ 

t  r* 

'  grr  4-  P*  ’{?*$•)  fp'f'*?) 

rt  / 
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W e  now  let 


with 


~  <f o  -^0  =  '}<»$,'?) 


H*-  (<*,?*.  f./?)  -  -£-^ 


'fc-i}0 


~*i !“»//.  rr  rrbl7  - 

Derfnrmpf)  _  “8 rai»  the  q0  integration  is 


performed  analytically  as 

£/'  Crt(fo.'?o.f.Q)  , 

bl>  COJ^(fa  -7e)  dric  «  ^  ^  *  /, 

where 


(27) 


c"  t  I  A . 


(28) 


^■l^^'^cVfr.n1 .  c|/-f,| -wr-r.w[, 

*  c  *<«  rJj)  [/>  c*5  /y  ^  c  ^  ^ 

4,  .  [\aa/{b,'^.i),,c '((. f, )'J I  .  c|  r.  f.\->tn(r-r.)  s.n  ,[»  c„ 

#  c  j4/>  ^  2"  f.)]\[6eej  i>(q- 1)  *  c sen 


T*  *  Tj^-Ai/I .  cstn  fi(f-fj\ 

•  f  St$  c'('/.f.)']h  b(rt.  ,)st„(Mtt) 

S. 

%  *  ■  eoVy  • ,}  •<»**«  t.)  | 


(30) 


Examination  of  J3  shows  that  it  is  singular  at  3  cot iftfa+l)  if 

-ts  f.  s  /  .  ^  is  singular  at  ft,.  fot*  if  -/«£,«/.  These 

singularities  must  be  isolated  before  the  £0  integration  can  be  completed. 
Then,  the  second  integral  in  Equation  (26)  is  evaluated  as 


bl  -fu  i  cos-fa.-!?.)  *df‘  ml  1 

■  •  i}jnr  ■  7*!jr  I 

y  z-/1  JJ*(€o.€.  •?)  d£o  (31) 


The  following  integrals  arc  evaluated  analytically: 

/  '?)<*£,  -  m  ^  jl?cas  f(^*0  *  C  Sin  bcos  4>(q*1) 

*-C  StnMf.f)  |  t(ljtf)  ♦  CJ//»  ^n/  #f)|  ♦  2c  S*n 

*  ^  n  * T* 0  -  c(f  \b *(q*l) **  *  b(q*l)  j 

-  MrOJ*2\(>'(9.»t+c,H'i)'\K *  c(r  o\  +  cU.fitJpfq.tT.t'k . *>r7.i)  | 

-6iq>t,Jn>2 1  c*(f  -  0'\^  *  c  (l-  f  «•  £o  J  (32) 

l !.((.  !.,) if..  J,u. ,) 

•  c*t»tu Mcost(^  t) 

-  t  J*(q  ()A+,lf>(t‘f)  -  c((»t)Jr%  .!)*•  €*((•* fl\  *  j 

'1 A  2  :(i  Vf  'Ac7'/ •  j.  c(ffU»  | 
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Also,  let 


J.  *  * 

-  HitZcjo) St 

^4  (fo.  £?)  m  Hi  (<•.'?*)  -  Hi  ( £at »  'jo  ) 


cos  V 

Then,  Equation  (26)  becomes 


(34) 


(35) 

(36) 


x- -£iZ t  ( '■% :)* 

Pm  /  /f  9  t 

or  /»./ L 

+  H,Ur.  f.rr)Ts({r.{.q)  *  ?,({>.*},)  Ir  (f.S.  ?)} 


0r\  TT^J” 


(37) 


which,  when  /^ar c  chosen,  is  in  (orm  (or  direct  numerical  evaluation. 

The  unknown  c.'t  are  contained  in  and  .  and,  be(orr 

solution  can  be  accomplished,  they  must  be  separated  so  that  Equation 
(37)  will  be  in  the  form  o(  Equation  (13).  I ,  is  then  treated  similarly 
to  f,  and  the  resulting  expression  is 


r, 

If  •  “  T~  Z  Z"  H,  (t>,  •?,(,  f) 

§rd>  /w/*/ 


‘#7  {  iv f)*st <*. *'?)**(&'*  fJ 

*  e  ().(  f )  7a  f  *  ^  (f  *7/  ] 

.  7-yjrjjjr .  '»«■*»/  •>*. 


whfff 


«  0  f  *  if*  <?#)  r  t,  t 

**  "  - *',t*f*) 

(J9) 

5,  )  w  -  £3L| 

*  /. 

HO) 

t  /#  4  %(€",.  f#  * 

' '* f-'  ’  .w. 

HU 

'/ ^  ^ 

Hi) 

f)  *  ,f,  )  •  S,  )(LllT 

"■V 

(4!) 

^  ■"  ^ 

/.f)  ♦  4rf&  y»)  ■ 

1*  4* 

(44) 

<*  -  /•  *  CAf?(f  f,) 

(45) 

?/,  *  f '  #  <«*  ^ 

(46) 

fv  *  T  *  *  f(im  *) 

(47) 

f 

«,(C.+.t.  f)  -  —*“*>//,  f,  < 

(41) 

<‘•141 

(49) 

1 1  1 


where 

A3  =  i>\il->{o\-c0s  <P  sgn(q-qe)\csLnt(f-l) 

+  bcos  |  [cst/i  <!>({+ 1)  *■  b  cos  J ?(?-/?*)] 

i 

A4  =.  f\un+\ \b Z(f] -q.)*+  c*f £+/)*]*  ^  ^  1 7  -  7„|  -costsgn(q-ije)[cstn PU-1) 

4-b  COS  -qe  )jj[c.As»  f'(f-t)  *■  b  cos  f  (q-qa) | 


(50) 


(51) 


_  b_ 

Ksfacf,  7)  -  ■j^^J~\csv>f(?-t)+l>COsf(r}-t ffc)  | 

-sqnCsinri—^Jn,  | [b^rf -q0^* c'fe- //*] 1 4 c/(-l) syn(sm p) |  ^52) 

_  b 

K*  (?,.£.?)  •—^T—‘6*\cS">f('f*‘f)  +  teosf{q-q0)\ 

t-Sfn(suif)  t-J*v  |  \b*(f]-f}o)9rC*(ft’1)*\l*  c( ?  *  1)  sfn(stnf)  j  ( 5  3  j 

*><?,./>).  >  -y  *  -?/]* 

S*'><f(?-q.)\c  Sm  Mf-  O+b  cos  f  (7 /},)] 

tf'?’?*)  +\c*(f'k/)e+bYf0 ■f)*]*’ 

Sin^(q-q0)\c Sin  t(f* t)  +  b  COS +(q- q, j] 

*(,(£■•})  «J  |f^  603  ft*?’*)  A  CScnf(f-f)\j'*\bcost(qf)*CSwt(ff)\ 

-\bcOs4>(qtf)t  csm  f(f-t)  Jn\bco$b(q*t)*  CSin+(f  f)\  *  /b  cost*  I 
-  [ic (t  {)U 2c(r-() - b(f , * »'* i(f- i/J1  *cO  !)\ 

+  ib  *c(t-f)Ut \[b'(f,t)‘ t c*(f‘/)']{ *  bff  * t) | 

+  b(q  I)hs.\b'(f  //‘fVf  O'j^eO  f)\  !)*],*•  k'l- f)|J  (55) 
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♦ 


*7  <£?)=•  ^bcos  ?ty-r)  +  csinil>fc+r)]j*o\bcos+(f-f)  *  cstn  ?  +f)  j 

-^bcos1>(q+l)*  csin  j>(f+f)\jn\b  COS  {(rf+l)  +  CSir>t(g+l)\  *  2b  COsf*  J 
+  ^^^\^2cO^)J^2cU^)-b(^l)J^\bz(^lf-c(t*lf\ -c(4+l)\  +2i 

4-o(£+0Jns2  \\b  (fftl)  4-C  U *0  J  +b(/f*f)\ *b(^-t)Jru  l[p*(rj-t)Z-bC*(f+1) 

-c(f+f)  | * c(f+f)£„ 2 1 [b’fy-n % c({+1) *+  b(i-n)  I  ]  (56) 

Here,  <6  and  K7  have  been  specialized  for  the  range  of  t  of  interest. 
For  convenience,  we  introduce  the  following  functions: 


<oo 


n- 


aV) 


-2, 


7a 


7o 


J  7o~7a 


Q  .  UA'ri)  -v.O-rJ)  _j-’  A  J_ 

f  7o~7a  (  J?#-tf# 


7a 


O,.  .  -ML±} 

L  7#  "  7/  )  7*  *  7a 


[d  in»  ini)  ~(i  "7  7# -ini)  ,  _  I 

i - FT, - *",WJ 


IIS 


f)  _  f(i  '  ifr  ~i£p)  -  (r~  TT&  ~T?*)  1 

- j - 7T^ - 2^*ta 


Mxc*°  for  \rjo,\z1 


MX7  *  O  for  |^J  2-  / 

*>»  -  */'-**;*  J 


II 

It 

“ 

r 

Not  ~ 

*r/- 

n03  - 

- 

3#/ 

*04  » 

*£T  ■ 

£ 

C* 

II 

-z(f 

Then  the  As  and  Is  can  be  expressed  as  follows; 

,  ,rt*  ~  j  ( t'fr)  ^oo  Jt 

+L(f +S(r-4;fr*4  J  J,  *-[(t-t,)Rco-2(i-z',){  r.,]jf 

*  t~is  {t-qf)Jr  ~4r  ^  ’srr^Af  ^7 


M<X.  - 


1-Vo, 


Mfe.  =  (’-&)* 
^2C  *  fad-j/,)* 


n 


07 


(f-nA)* 


**t7  * 
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JL  ^  (l~£p)  @OfJl  ~  a  j  ( t  ~  $ p)  Rot  J~2 

P*f  fi*/  Gr  P*f  i 

+ 1 \U-Sp)Rot  rae]J3  +  \(i-z,)R0,  -2(f-ft){r07]  j4 

+  (l~£p)(  t-tfo)  foZrj  Te?  27 

Z  f  \(f-fr)Xo2Jz 

P*t  fi=t  or  p=t  l 

~4p)  Rot  ~  2(}-<f/>  )*  To&  J  Ts  +\(.1-%p)Ro2.~2('I~$p)^  T~o7 

+  (1-fr)(l-f?o‘)ir.*Ts}  “  ToiJ,  -±  T07fo7 


Z  Z  (f-ltfat-fr)  Q*S  J>  ~jk  Z  \(f-*r)  *"  J* 

p*i  js~f  or  2*f  >■ 

+  [('/-/>;  *03  -20-&)*  To^Js  t[('-fr)  Zo3-2(t-ti?r07]  J4 

+  0-  *A  (’-V)  io'J,  j  -£  T„  J,  -  ±  r0?  Jr 

Z  X  ('•&)*  (*’M  Go*  *  -rl  \(i-?p)*o*K 

psf  flzf  °r  p-  1  L 

+  [(1-ip)  R04  *-2(f-fp)i  T0i}T3+[(l-{p)Ro4  -2(1  -ip)'  r07]  J+ 

+  (l~£p)  (I'jo)  ^0  ~  y~  TOi  J*  *  T07  ^7 


Z  Z  (  l-'ffi  )*  ( 1~fp)  Qoo  <?/  -  z~  Z  f  (fmtp)  ^ °° 

P*t  fl*r  r  p*t  l 


+  [( 1~  ip  )  Rco  +  2  ( /- €p)*  +  \(l'ir  )Roo -2(t~fp  )*  ^7] 


j; 


+  ( t-  )  ( 1-P~)  Tf  ~l  -  - —  TtL  J*  4-  .  Tj ,  Jr 


a"  -  -£:Z  j (>-&)**, s* 

arA'p*tfl*t  8r  £rf  t 

r  [( 1  ~  Zp)  Roi  -  2(l-fZ)i  7>J  T3  +  [( 1-$p)Ro,  -<?<?'<£>)*  r,7 ]  I4 

*  (t~Rp )(  )  ?<>  Jj~  J  —  Ttl)  Tf7  u~7 

f*  AS  y> 

*«  -  -j-z 

P*  /  fi*1  or  /»*  /  L 

+  [< t-fr )  *02  +2(1-*?)*  Ttc J  Jj  +  )  /p^  -  2(t-&)i  r„]  J4 

+  (f’&X'-j,*)?:*}  "  “  T,uTb  t-~  7/7  J" r 

or'd/  Pit  fi*f  9r  p7f  (■ 

+  [0-  fp)  *03  -  2  ('-{?)*  T«]  J3  V-  [( 1-Zp)R05 T,j]  J4 

+0-4p)(l-no)no^s]  -£rT”Jr 

At4'  \(l'*r)*o*Tz 

P=tj*/  8rpTt  <- 
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